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Abstract 


The  Air  Force  Weather  Agency  (AFWA)  has  two  operational  algorithms  that  derive 
surface  temperatures  from  microwave  observations  taken  by  the  Special  Sensor 
Microwave  Imager  (SSM/I)  which  rides  aboard  space  platforms  of  the  Defense 
Meteorological  Satellite  Program  (DMSP).  The  algorithm  called  Temperature  from 
Satellite  Microwave  Imager  (TS)  is  used  to  forecast  surface  temperatures  and  analyze 
global  cloud  coverage.  The  second  algorithm  is  fittingly  called  CalibrationA^alidation 
(CV),  as  it  was  the  algorithm  used  to  calibrate  and  validate  the  first  SSM/I  back  in  1987. 
Multiple  linear  regression  defined  the  algorithms  from  empirically  gathered  brightness 
temperatures  and  simultaneous  surface  temperatures.  The  key  questions  at  hand  are  how 
much  data  do  these  algorithms  produce  and  how  accurate  is  it. 

The  current  thesis  answered  the  above  questions  with  a  multi-seasonal  comparative 
study  over  four  locations.  The  study  matched  algorithm  outputs  to  conventional  weather 
station  temperature  readings  to  produce  data  pairs.  Over  13,300  data  pairs  were 
generated  from  the  1996  summer  and  fall  and  1996-1997  winter  for  the  Continental 
United  States  (CONUS),  Bosnia,  Korean  Peninsula,  and  Saudi  Arabia. 

The  results  show  TS  produced  on  average  7  %  more  surface  temperatures  than  CV. 
CV  met  AFWA’s  accuracy  criteria  16  %  more  often  than  TS.  However,  CV  was  only  1.0 
degree  Celsius  more  accurate  than  TS  on  average.  The  study  also  generated  bias  tables 
for  all  locations  and  seasons. 
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EVALUATION  OF  SATELLITE  MICROWAVE  DERIVED 


SURFACE  TEMPERATURE  ALGORITHMS 
FOR  THE  PERIOD  AUGUST  1996  TO  FEBRUARY  1997 


I.  Introduction 


Chapter  Overview 

After  a  general  introduction  to  the  problem  of  inadequate  amounts  of  global  surface 
temperatures,  the  chapter  introduces  the  concept  of  remote  sensing,  specifically 
microwave  remote  sensing.  The  background  section  outlines  the  space  platform, 
instrument,  and  algorithms  of  concern  in  this  study.  After  the  specific  research  problem  is 
defined  and  some  assumptions  made,  the  scope  and  general  approach  of  the  research  are 
explained.  Finally,  this  chapter  provides  the  general  results  of  the  research. 

Introduction 

One  special  attribute  of  the  human  species  is  the  ability  to  extend  our  natural  five 
senses,  especially  sight,  with  tools  of  our  own  creation.  By  creating  sensors  that  detect  or 
see  beyond  the  visible  part  of  the  electromagnetic  spectrum,  we  learn  more  about  our 
environment.  One  class  of  these  sensors  resides  on  earth-orbiting  satellites;  they  detect  a 
totally  different  portion  of  the  electromagnetic  spectrum  -  microwaves.  These  microwave 
observations  determine  many  atmospheric  parameters,  including  surface  temperature.  As 
a  result,  humans  extend  their  natural  senses  to  see  global  surface  temperatures  from  space. 
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Background 

Even  in  modern  times,  the  ability  to  accurately  measure  global  surface  temperatures  is 
a  problem.  The  conventional  weather  station  network  does  not  cover  remote  regions. 

The  data  sparseness  affects  atmospheric  models  that  require  a  fine  data  analysis  to 
accurately  predict  the  future  state  of  the  atmosphere.  As  pointed  out  by  McClatchey  and 
Greenwood  (1996),  the  timeliness,  specificity,  and  responsiveness  requirements  of  the 
military  for  weather  information  greatly  exceed  those  of  the  civilian  world.  Remote 
sensing  from  space-borne  platforms  is  one  answer  to  the  data  sparseness  problem. 

One  method  to  improve  the  analysis  of  global  surface  or  land  temperatures  is  through 
the  remote  observation  of  earth’s  microwave  emissions  sensed  from  satellites  -  a  process 
referred  to  as  Passive  Microwave  Radiometry  (PMR).  Space-bome  PMR  has  been  in 
existence  since  the  early  1960’s  (Ulaby,  1981).  The  primary  advantage  of  using  the 
microwave  portion  of  the  electromagnetic  (EM)  spectrum  (versus  thermal  or  visible)  is  the 
ability  to  see  through  clouds. 

There  are  many  types  of  remote  sensing  satellites.  One  such  system  is  the  Defense 
Meteorological  Satellite  Program  (DMSP).  The  type  of  microwave  sensor  aboard  the 
DMSP  satellites,  called  Special  Sensor  Microwave/Imager  (SSM/I),  has  been  operational 
since  1987.  In  April  1997,  the  US  Air  Force  launched  the  latest  DMSP  satellite,  F14, 
from  Vandenberg  AFB,  California  (Cooper,  1997). 

Air  Force  Weather  Agency  (AFWA)  has  two  operational  algorithms  for  calculating 
surface  temperatures  derived  from  microwave  emissions  sensed  aboard  DMSP  space 
platforms.  The  first  algorithm,  called  Temperature  from  Sensor  Microwave  Imager 
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(TMPSMI  or  TS),  provides  input  for  the  global  surface  temperature  model  and  the  global 
cloud  coverage  model  (Kopp,  1995).  The  second  algorithm,  fittingly  called 
CalibrationA^alidation  (CV),  was  the  algorithm  used  to  calibrate  and  validate  the  first 
SSM/I  in  1987. 

Problem  and  Assumptions 

AFWA’s  atmospheric  models  require  an  accurate  global  analysis  of  surface 
temperatures  in  order  to  correctly  specify  and  forecast  global  cloud  coverage.  In  addition, 
military  forecasters  need  near  real  time  surface  temperatures  from  data  sparse  or  denied 
territory.  A  current  temperature  analysis  aids  in  model  verification  and  temperature 
forecasts.  Among  many  other  military  uses,  these  surface  temperature  applications  are 
critical  for  the  use  of  precision  guided  munitions  -  key  tools  in  the  application  of  today’s 
US  military  power. 

The  problem  addressed  in  this  paper  is  the  performance  of  the  TS  and  CV  algorithms. 
Which  algorithm  produces  more  temperatures  and  which  one  is  more  accurate?  No 
research  to  date  has  extensively  compared  the  two  algorithms. 

There  are  two  assumptions  from  the  outset.  First,  the  weather  station  temperature 
observations  are  truth  data.  There  was  no  attempt  made  to  analyze  the  accuracy  of  the 
station  data.  Secondly,  the  algorithms  contain  correct  computer  code.  This  study  does 
not  aim  to  improve  the  algorithms,  rather  just  compare  their  outputs. 
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Research  Scope  and  General  Approach 

DMSP  SSM/I  data  from  the  F13  satellite  was  the  only  microwave  data  source  in  the 
study.  An  initial  study  occurred  in  the  fall  of  1996  but  portions  needed  reaccomplishment 
due  to  data  problems.  Furthermore,  the  comparison  needs  spatial  and  temporal  expansion 
(Comoglio,  1997). 

SSM/I  and  weather  station  data  were  gathered  for  summer  1996,  fall  1996,  and  winter 
(Jan-Feb)  97.  The  study  covers  the  land  regions  of  the  CONUS  and  other  regions  of 
interest  to  the  US  military:  Korea,  Saudi  Arabia,  and  Bosnia. 

After  reformatting  the  raw  SSM/I  data,  both  algorithms  were  run  for  the  areas  of 
interest.  The  production  counts  were  then  obtained.  The  raw  station  data  was  also 
reformatted.  Then,  a  program  matched  the  SSM/I  and  weather  station  data  by  location 
and  time  to  create  data  pairs.  Finally,  descriptive  statistics  were  generated  and  analyzed 
on  production,  tendency,  and  accuracy  for  all  regions  and  seasons. 

This  study  does  not  present  information  as  to  how  DMSP  satellites  are  tested  for 
sensor  accuracy.  The  process  to  ensure  sensor  readings  are  correct  (e.g.  calibration)  is  left 
for  others  to  evaluate.  However,  error  possibilities  will  be  presented. 

Results 

The  study  created  a  combined  total  of  over  13,300  data  pairs  for  twelve  data  sets:  the 
four  locations  mentioned  above  for  summer,  fall  and  winter  seasons.  Overall,  TS 
produced  7  %  more  data  than  CV.  For  most  locations  both  algorithms  exhibited  a  cold 


4 


bias:  the  algorithm  temperature  had  the  tendency  to  be  colder  than  the  actual  temperature. 
However,  CV  had  a  warm  bias  for  all  three  seasons  in  Saudi  Arabia.  CV  met  the  accuracy 
criteria  of  <  8  degrees  Kelvin  (K)  Root  Mean  Square  Error  (RMSE)  16  %  more  often 
than  TS.  However,  CV  was  on  average  1.0  K  RMSE  more  accurate  than  TS. 

Summary 

The  main  purpose  of  this  research  is  to  provide  the  production  and  accuracy 
characteristics  of  the  TS  and  CV  algorithms  for  AFWA.  A  comparison  study  showed  how 
the  algorithms  performed  in  different  geographic  locations  and  at  different  times  of  the 
year.  The  algorithm  outputs  were  compared  to  each  other  as  well  as  to  the  conventional 
weather  station  surface  temperatures.  Overall,  TS  produced  more  output  than  CV,  but 
CV  was  more  accurate. 
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II.  Backoround 


Chapter  Overview 

After  an  introduction  to  the  physics  of  electromagnetic  radiation,  this  chapter  provides 
a  short  history  of  microwave  remote  sensing.  The  satellite  platform  and  its  microwave 
sensor  particular  to  this  study  are  then  discussed.  Special  emphasis  is  placed  on  what  the 
sensor  sees  or  detects.  In  turn,  what  the  microwave  sensor  provides  as  output  is  discussed 
in  the  context  of  antenna  theory.  A  radiative  transfer  equation  describing  radiation 
received  at  the  satellite  is  then  outlined  in  detail.  Some  potential  solutions  to  such  an 
equation  are  reviewed,  including  the  solution  method  of  choice:  Statistical  Inversion 
Method.  Such  a  solution  results  in  the  two  algorithms  under  scmtiny  which  are  compared 
and  contrasted.  After  discussing  applications  of  the  microwave  derived  surface 
temperatures,  the  chapter  presents  sources  of  error  in  retrieving  surface  temperatures  from 
microwave  measurements.  Finally,  the  chapter  ends  with  a  current  literature  review  of  the 
determination  of  land  surface  temperature  via  remote  sensing. 

Introduction 

The  algorithms  under  consideration  use  microwave  radiation  sensed  aboard  satellites  to 
calculate  “skin”  temperatures.  “Skin”  refers  to  the  very  top  layer  of  the  bare  soil, 
vegetated  soil,  and/or  snow  cover.  This  study  compares  the  skin  temperatures  to 
conventional  weather  station  surface  temperatures.  The  standard  height  of  a  thermometer 
at  a  weather  station  is  1.5  meters  above  the  ground  (Mcllveen,  1986).  For  this  thesis. 
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both  skin  and  weather  station  temperature  will  be  referred  to  as  surface  temperature,  but 
conceptualized  as  the  temperature  of  the  air  near  the  ground.  Before  describing  the 
algorithms  that  calculate  the  land  surface  temperature,  much  background  information  will 
be  covered  to  fully  understand  how  the  algorithms  function. 


Electromagnetic  (EM)  Radiation 

All  objects  not  at  absolute  zero  (0  K)  emit  energy  across  a  continuous  range  of 
frequencies  or  wavelengths.  The  entire  range  of  frequencies  comprises  the 
electromagnetic  spectrum  (Rees,  1990). 
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Figure  1.  Wavelength  Distribution  of  Earth  (T=300  K)  with  Highlighted  Microwave 
Region.  B(X.)  =  Brightness  (Watt  m'^  cm'*)  and  X  =  wavelength  (cm). 
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Each  object  has  a  particular  distribution  of  wavelengths  across  the  EM  spectrum  that  is 
a  function  of  the  object’s  temperature.  The  focus  here  is  the  earth- atmosphere  interface 
where  the  surface  temperature  is  measured.  The  wavelength  distribution  for  such  an 
interface  is  shown  in  Figure  1.  The  figure  presents  the  amount  of  the  radiation  emitted, 
termed  brightness,  plotted  against  wavelength  using  a  mean  surface  temperature  of  300  K. 

The  peak  of  the  curve  in  Figure  1  shows  the  earth-atmosphere  interface  primarily 
transmits  in  the  infrared  portion  of  the  electromagnetic  spectrum  (3  |4.m  to  15  p.m).  Most 
weather  satellite  measurements  observe  the  infrared  (IR)  and  visible  (VIS)  portion  of  the 
EM  spectrum,  the  most  common  being  the  satellite  imagery  used  everyday  in  TV 
meteorology  broadcasting.  However  when  sensing  surface  temperatures  from  space,  the 
IR  or  VIS  portion  of  the  EM  spectrum  is  fraught  with  difficulties  as  other  radiation 
sources  interfere  with  the  emission  from  the  surface.  Thus,  scientists  look  to  other 
portions  of  the  EM  spectmm.  Figure  1  also  shows  a  significant  but  smaller  portion  of  the 
earth’s  emitted  energy  is  in  the  microwave  region  (vertical  markers).  The  microwave 
region  is  generally  considered  to  span  the  range  from  0.1  cm  to  10  cm,  or  in  terms  of 
frequency,  from  3  GHz  to  300  GHz  (Janssen,  1993). 

The  microwave  region  looked  promising  for  several  reasons.  Ulaby  (1981)  gives  four 
main  reasons  for  the  use  of  microwaves  in  remote  sensing.  First,  microwaves  can 
penetrate  clouds  (unlike  IR).  Second,  microwaves  are  independent  of  sun’s  illumination 
so  that  the  sensor  can  measure  during  the  nighttime  (unlike  VIS  but  similar  to  IR).  Third, 
microwaves  are  able  to  penetrate  into  vegetation  (such  as  dense  jungle)  and  even  soil. 
Finally,  microwaves  are  simply  different  from  VIS  and  IR .  Data  gathered  in  all  three 
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regions  of  the  EM  spectrum  would  complement  each  other  and  yield  more  information  in 
the  aggregate. 

This  study  compares  the  surface  temperatures  observed  at  weather  stations  to  satellite 
derived  temperature  readings  of  the  same  portion  of  the  atmosphere.  What  follows  is  a 
contrast  of  the  processes  and  instruments  used  for  those  two  readings. 

The  standard  bulb  thermometer  at  a  weather  station  measures  the  surface  temperature 
through  conduction,  i.e.,  contact  with  the  air.  The  molecular  motions  (kinetic  energy)  of 
the  air  molecules  are  transferred  to  the  liquid  molecules  in  the  thermometer.  In  contrast, 
the  geometric  (orientation)  and  emission  properties  of  the  earth’s  surface  and  intervening 
volume  (atmosphere)  determine  the  microwave  temperature  readings  from  an  earthward 
looking  satellite  (Ulaby,  1981).  The  latter  process  is  radiative  transfer  of  energy  and  the 
former  process  is  conductive  transfer  of  energy.  For  space-based  observations,  radiation 
and  the  effects  during  the  radiation’s  propagation  through  the  atmosphere  determine  the 
observed  surface  temperature.  As  such,  thermometers  are  direct  sensors  while  space 
observations  are  remote  sensors.  “Radiometer”  is  the  name  given  to  a  remote  sensing 
instrument  that  measures  EM  radiation.  A  radiometer  is  a  receiver  (passive  system)  as 
opposed  to  radars  that  send  as  well  as  receive  EM  energy  (active  systems).  The  remote 
sensing  instrument  of  interest  is  referred  to  as  a  passive  microwave  radiometer  and  is  a 
very  sensitive  receiver  capable  of  measuring  low  levels  of  microwave  radiation  (Ulaby, 
1981).  Figure  1  shows  the  relatively  low  amounts  of  microwave  radiation  emitted  from 
the  earth-atmosphere  interface. 
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With  the  goal  of  understanding  how  passive  microwave  radiometers  detect  surface 
radiation,  the  emitted  radiation  from  the  earth-atmosphere  interface  must  be  described. 

The  flow  of  radiation  from  a  source  to  a  receiver  is  called  radiative  transfer.  An  equation 
describes  the  transfer  process  -  how  the  radiation  is  affected  from  source  to  receiver. 
Three  main  effects  need  explanation:  ground,  atmosphere  (medium),  and  radiometer 
(receiver).  The  radiative  transfer  equation  can  have  many  forms.  The  equation  developed 
here  describes  microwave  energy  flowing  from  the  earth  and  sky  through  the  atmosphere 
until  it  reaches  the  satellite  receiver. 

Assuming  the  earth  had  no  atmosphere  and  was  a  perfect  absorber  and  emitter  of 
radiation  (referred  to  as  a  blackbody),  Planck’s  Law  describes  the  amount  of  emitted 
radiation.  Planck’s  law  states  the  surface  brighmess  is  the  flow  of  energy  across  a  unit 
area,  per  unit  time,  per  unit  frequency  interval,  from  in  all  directions,  and  is  given  by 
(Ulaby,  1986): 


^2  ehv/kT_j 

where  Bv  =  Surface  brightness  or  radiance  (W  m'^  sr‘Hz'') 
V  =  Frequency  (Hz) 

T  =  Physical  temperature  (K) 
h  =  Planck’s  constant  (6.6256x10'^  J  s) 
k  =  Boltzman’s  constant  (1.3805x10'“  J  K'*) 
c  =  Speed  of  light  (2.9979x10*  m  s  ') 


(1) 


Note  Equation  (1)  is  for  a  particular  frequency.  The  radiometer  of  concern  in  this 
study  uses  four  frequencies  (more  on  this  later). 
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Since  the  focus  here  is  the  microwave  region,  the  Rayleigh- Jeans  approximation  can  be 
used.  An  exponential  Taylor  series  is  applied  to  estimate  the  e*  term  in  Equation  (1) 
(Spiegel,  1993  &  Ulaby,  1981) 

x^  x^ 

e*  =l  +  x  +  — -f-— (2) 

2!  3! 


With  X  =  hv/kT  and  noting  for  microwaves  hv/kT  «  1  for  Earth-like  temperatures, 
Equation  (2)  simplifies  to  =  1+x.  In  turn.  Equation  (1)  simplifies  to: 


BJT)  = 


2v^kT 


(3) 


or  alternatively  since  \  =  cfk.  Equation  (3)  can  be  changed  to  the  wavelength  form: 


2kT 

A' 


(4) 


The  important  outcome  of  the  Rayleigh- Jeans  approximation  is  the  linear  relationship 
between  brightness  and  temperature.  Figure  2  shows  the  linear  approximation  to  the 
original  Planck  function  is  quite  good  for  the  microwave  region.  In  fact,  the  Rayleigh- 
Jeans  approximation  has  only  a  1  %  deviation  from  Planck’s  law  for  earth  applications  if  X 
>  2.57  mm  or  v  <  1 17  GHz  (Ulaby,  1981).  The  radiometer  used  in  this  study,  called  the 
Special  Sensor  Microwave  Imager  (SSM/I),  functions  within  these  constraints. 
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Non-Blackbody  and  Land  Effects 

A  blackbody  is  an  idealized  material  that  absorbs  all  the  incident  radiation  at  all 
frequencies,  reflecting  none.  Another  aspect  of  a  blackbody  is  it  acts  as  a  perfect  emitter. 
It  is  assumed  the  blackbody  is  in  thermodynamic  equilibrium  with  its  environment:  energy 
is  absorbed  and  emitted  at  the  same  rate  resulting  is  a  constant  temperature. 
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Figure  2.  Rayleigh- Jeans  Approximation 
(Brj  &  dashed  line  on  curve)  to  the  Planck  Function. 

Real  materials,  non-blackbodies  such  as  the  earth’s  land  surfaces  and  atmosphere,  emit 
less  than  and  do  not  absorb  as  much  as  a  blackbody.  The  physical  property  called 
emissivity  accounts  for  this  departure.  Emissivity,  e,  is  defined  as  the  ratio  of  brightness 
of  a  non-blackbody  (B)  to  the  blackbody  brightness  (Bbb)  or  ratio  of  the  brightness 
temperature  (Tb)  to  the  physical  temperature  (T)  for  microwave  radiation  (Ulaby,  1981): 
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(5) 


In  other  words,  emissivity  is  a  measure  of  a  surface’s  ability  to  radiate  energy  and 
depends  on  the  molecular  structure  of  the  surface  (Conway,  1997).  Tb  varies  from  zero 
up  to  a  maximum  value  equal  to  the  physical  temperature  of  the  scene.  Representing 
radiation  from  a  scene,  Tb  is  the  end  product  of  a  radiometer  measurement.  However,  T, 
the  physical  temperature,  is  the  goal  of  algorithms  in  this  study.  Exactly  how  the  physical 
temperature  is  retrieved  from  the  brightness  temperature  will  be  explained  in  later  sections. 
Equation  (5)  assumes  the  material  is  homogeneous  and  of  uniform  temperature.  Both 
assumptions  may  be  an  oversimplification  for  sensing  temperatures  near  the  ground.  In 
contrast,  the  assumption  of  homogeneity  is  more  valid  for  an  ocean  surface.  As 
McFarland  (1991)  points  out,  land  emissivities  are  more  variable  both  in  space  and  time. 
Spatial  variances  arise  from  vegetation,  topography,  urban  features,  bodies  of  water,  etc. 
Examples  of  temporal  changes  are  caused  by  vegetation  changes  both  natural  and  artificial 
(agricultural),  soil  moisture  changes,  and  snow  cover  morphology.  For  both  land  and 
ocean  surfaces,  another  emissivity  complication  is  surface  roughness,  which  generally 
reduces  the  radiance  sensed  by  the  radiometer  by  diffusion  (Ulaby,  1993). 

Another  surface  property  at  the  earth-atmosphere  interface  is  reflectivity,  r.  Some 
energy  reaching  the  radiometer  is  due  to  reflection.  Reflectivity  is  closely  related  to 
emissivity.  Assuming  no  transmission  through  the  body,  non-blackbodies  reflect  the 
incident  energy  if  it  is  not  absorbed,  as  shown  by  Equation  (6). 


13 


1  =  e+r 


(6) 


Reflectivity  at  a  surface  is  also  referred  to  as  siuface  scattering.  The  equation  for 
emissivity  (Eq.  5)  implies  if  the  brightness  temperature  measured  by  the  satellite 
radiometer  and  the  surface  emissivity  are  known  then  the  physical  temperature  of  the 
scene  can  be  determined  assuming  no  atmospheric  influence  on  the  radiation  as  it 
propagates  between  surface  and  receiver.  This  is  one  approach  in  determining  surface 
temperatures  from  space.  Unfortunately,  several  complexities  arise  which  are  explained  in 
later  sections.  Another  aspect  of  EM  radiation  is  polarization. 

Electromagnetic  radiation  travels  in  waves.  As  the  name  implies,  two  fields  define 
radiation:  electric  and  magnetic.  A  combined  wave  simultaneously  represents  each  of 
these  fields.  These  two  waves  are  always  perpendicular  to  each  other  and  to  the  direction 
of  propagation.  Polarization  is  the  orientation  of  these  fields  (Rees,  1990)  with  respect  to 
a  coordinate  system.  For  example,  a  sensor  could  be  sensitive  to  the  vertical  or  horizontal 
orientation  of  the  EM  wave. 

With  the  basic  EM  knowledge  and  terminology  known,  the  next  section,  paraphrased 
from  Ulaby  (1981),  provides  a  context  of  the  present  state  of  passive  satellite  microwave 
radiometry  (PMR). 

Remote  Sensing  History 

Using  frequencies  near  the  microwave  portion  of  the  EM  spectrum,  the  first  remote 
sensing  experiment  occurred  more  than  1 10  years  ago  in  1886.  In  the  1920’s,  the  Naval 
Research  Lab  (NRL)  pioneered  the  use  of  radars  to  detect  ships,  aircraft,  and  the  height  of 
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the  ionosphere.  The  NRL  is  mentioned  because  they  served  as  technical  consultant  to  the 
Department  of  Defense’s  (DoD)  atmospheric  remote  sensing  system  under  investigation  in 
this  paper.  Airborne  detection  systems  using  microwaves  became  a  reality  in  the  early 
1940’s.  As  opposed  to  the  active  systems  just  mentioned,  PMR  first  developed  in  the 
1930’s  and  1940’s  to  measure  radiation  from  space.  Terrestrial  microwave  radiometry 
(looking  down)  had  its  beginning  in  the  late  1950’ s.  Synthetic  aperture  radar  began  in  the 
1950’s.  This  was  a  major  improvement  in  resolution  that  paved  the  way  for  space-borne 
imaging  radars  (Ulaby,  1981). 

NASA  deployed  its  first  space-borne  passive  microwave  radiometer  over  Venus  with 
the  Mariner  2  spacecraft  in  December  1962.  The  USSR  flew  the  first  earth-viewing 
atmospheric  passive  microwave  radiometer  in  1968  aboard  the  Cosmos  243  satellite.  The 
US  National  Oceanographic  and  Atmospheric  Administration  (NOAA)  followed  in  1972 
with  an  atmospheric  passive  microwave  radiometer  aboard  the  Nimbus  5  satellite.  The  US 
Defense  Meteorological  Satellite  Program  (DMSP)  deployed  its  first  atmospheric  passive 
microwave  radiometer  in  1978  to  sense  atmospheric  profiles  of  temperature  (Ulaby, 

1981).  The  Special  Sensor  Microwave  Imager  (SSM/I)  instrument  on  board  the  DMSP 
satellite  F8  was  first  sent  aloft  in  1987,  a  full  25  years  after  Mariner  2’s  pioneering  journey 
(Hesser,  1995). 

Atmospheric  PMR  has  concentrated  on  soundings  and  ocean  parameters  such  as  wind 
speed  over  water  and  sea  ice.  Researchers  investigated  some  land  parameters  such  as  soil 
moisture,  but  surface  temperature  has  been  neglected  in  microwave  remote  sensing 
(McFarland,  1991).  Surface  temperatures  are  important  to  all  the  military  services  for  the 
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following  reasons.  First,  most  atmospheric  numerical  models  use  global  surface 
temperature  as  a  key  input.  Second,  near  real  time  surface  temperatures  over  enemy  held 
or  data-sparse  areas  are  needed  for  target  weather  analysis  and  forecasts. 

The  pioneering  instmments  mentioned  above  were  rather  simple  sensors  using  few 
channels  with  limited  resolution.  As  of  this  writing,  the  DMSP  SSM/I  still  represents 
some  of  the  most  advanced  technology  in  space-bome  PMR.  Only  DMSP  SSM/I  data 
were  used  in  this  study.  The  next  section  details  information  about  this  system. 


DMSP  and  the  SSMI 

First  launched  in  1970,  DMSP  Block  V  satellites  represented  state  of  the  art 
technology  with  such  advancements  as  0.62  km  visual  imagery  resolution  -  the  highest 
resolution  of  any  weather  satellite  then  and  even  today  almost  30  years  later  (Fett  et  al, 
1997).  Currently  operational,  DMSP  satellites  F13  and  F14  are  in  a  near  circular  polar 
orbit  at  a  mean  altitude  of  850  km.  The  spacecraft  is  about  the  size  of  an  average  car  (3.7 
m  long  by  1.2  m  high)  and  has  a  design  life  span  of  two  years.  The  satellite  makes  one 
orbit  in  102  minutes  (Deuel,  1996). 

Figure  3  shows  both  the  ascending  and  descending  paths  of  the  SSM/I  swaths  over  a 
24  hour  period  (Janssen,  1993).  The  swath  width  is  1394  km.  Note  the  data  gaps  in  the 
lower  latitudes.  The  SSM/I  is  a  seven  channel,  four  frequency,  linearly  polarized,  passive 
microwave  radiometer.  The  seven  channels  are  sensitive  to  vertically  and  horizontally 
polarized  radiation  at  channel-center  frequencies  of  19.3  GHz,  22.2  GHz,  37  GHz,  and 
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Figure  3.  SSM/I  Global  Coverage 


85.5  GHz  frequencies,  excluding  the  horizontal  polarization  at  22.2  GHz.  The  sensor 
measures  land,  atmospheric,  and  ocean  surface  brightness  temperatiues.  In  general, 
imager  radiometers  use  frequencies  in  the  atmospheric  windows  to  monitor  surface 
features  such  as  the  ones  just  mentioned.  In  contrast,  sounder  radiometers  use  frequencies 
in  the  less  transparent  portion  of  the  EM  spectrum  to  capture  profiles  of  atmospheric 
parameters. 

The  four  frequencies  of  the  SSM/I  are  indicated  with  an  “S”  in  Figure  4.  The 
corresponding  wavelengths  for  19.3  GHz,  22.2  GHz,  37  GHz,  and  85.5  GHz  are  1.55  cm, 
1.35  cm,  0.81  cm,  and  0.35  cm,  respectively.  Figure  4  highlights  the  sampling  intent  of 
the  SSM/I.  The  19.3  GHz  channel  was  chosen  to  sample  the  lowest  layers  of  the 
atmosphere  as  the  atmosphere  is  quite  transparent  at  that  frequency.  The  22.2  GHz 
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(cm) 


Figure  4.  Percentage  Transmission  Through  Clear  Atmosphere  in  Vertical  Direction 
with  SSM/I  Frequencies  (S)  (adapted  from  Ulaby,  1986) 


frequency  takes  advantage  of  a  water  vapor  absorption  band.  The  37  and  85.5  GHz 
frequencies  were  chosen  to  take  advantage  of  their  respective  windows.  However,  the 
addition  of  the  85  GHz  channel  is  a  unique  advancement  over  previous  microwave 
imagers.  With  its  high  sensitivity,  the  85  GHz  channel  is  primarily  used  for  measuring 
rainfall  (Janssen,  1993). 

Located  on  top  of  the  satellite  (see  Figure  5),  the  instrument  consists  of  a  45  °  tilted 
parabolic  reflector  (24  x  26  inches)  which  receives  the  scene  radiation  and  feeds  it  into  a 
broadband,  seven  port  horn  antenna  which  then  feeds  into  the  radiometer.  This  whole 
assembly  is  oriented  normal  to  the  earth’s  surface.  Due  to  the  curvature  of  the  earth’s 
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X 


Figures.  SSM/I  Scan  Geometry 


surface,  the  viewing  angle  (0)  remains  at  a  constant  53.1  degrees.  In  order  to  scan  the 
scene  on  the  earth,  the  assembly  rotates  once  every  1.9  seconds  to  make  a  conical  scan 
angle  of  102.4  degrees  or  1394  km  on  the  surface  of  the  earth.  Figure  5  also  shows  how 
each  frequency  is  sampled.  On  scan  A,  the  19.3  GHz  horizontal  and  vertical  (H&V),22.2 
GHz  V,  and  37  GHz  H&V  channels  are  each  sampled  64  times  across  the  swath  width. 
Each  sample  is  1.6  degrees  in  width.  Simultaneously,  the  85.5  GHz  H&V  channels  are 
each  sampled  over  128  equal,  0.8-degree  intervals.  On  scan  B,  the  SSM/I  samples  only 
the  two  85.5  channels,  but  at  the  same  equal,  0.8-degree  intervals  128  times  across  the 
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scan.  Macelloni  et  al  (1994)  states  the  ground  resolution  of  the  19.3  GHz,  22.2  GHz,  37 
GHz,  and  85.5  GHz  frequencies  are  43, 40, 29,  and  13  km  respectively.  These  resolutions 
also  represent  pixel  sizes.  The  13  km  resolution  of  the  85  GHz  channel  has  one  of  the 
highest  resolutions  ever  for  a  microwave  radiometer  (Janssen,  1993).  Although  the 
frequencies  were  chosen  for  the  determination  of  various  atmospheric  parameters,  all  four 
frequencies  are  used  in  the  surface  temperature  algorithms. 

The  basic  functions  of  the  SSM/I  instrament  have  been  reviewed.  The  next  section 
explains  what  the  SSM/I  actually  senses  in  the  context  of  atmospheric  effects  and  the 
resulting  radiative  transfer  equation. 


Atmospheric  Effects 

What  the  instrament  detects  is  a  complex  combination  of  several  radiation  sources  and 
radiation-matter  interactions.  A  radiometer  looking  at  the  earth  senses  a  brightness,  the 
origins  of  which  are  diverse.  Shown  in  Figure  6,  the  energy  incident  on  the  sensor  is 
primarily  from  four  sources:  Land  emissions  (1),  land  reflections  of  downward 
atmospheric  entissions  (2),  land  reflections  of  downward  cosmic  emissions  (3),  and 
upward  atmospheric  emissions  (4)  (Ulaby,  1981). 

With  Figure  6  in  mind,  the  scene  that  the  radiometer  views  is  actually  a  cylinder 
pointed  in  a  particular  direction.  The  direction  is  quantified  by  an  azimuth  angle  ((]))  and 
elevation  angle  (0).  Sources  (1)  and  (3)  are  relatively  simple  to  account  for  in  that  the 
atmosphere  reduces  the  radiation  along  the  way  to  the  radiometer.  Sources  (2)  and  (4), 
atmospheric  emissions,  are  more  complex  because  each  layer  in  the  atmosphere  emits 
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radiation  which  is  simultaneously  reduced  by  absorption  before  it  reaches  the  receiver. 
Before  these  sources  and  effects  are  examined,  some  general  radiative  transfer  process 
must  be  described. 


Satellite  Radiometer 


Figure  6.  Satellite  Radiometer  Sources  of  Radiation. 

There  are  two  processes  that  occur  between  radiation  and  the  medium  (in  this  case  the 
atmosphere)  in  which  the  radiation  propagates:  emission  and  extinction. 

Extinction  of  radiation  by  the  medium  involves  the  two  effects  of  scattering  and 
absorption.  Scattering  here  refers  to  radiation  diverted  out  of  the  radiometer's  path.  In 
contrast,  emission  scattering  means  scattering  into  the  radiometers  view  path.  Absorption 
means  the  radiation  is  transformed  to  other  forms  of  energy,  such  as  heat  (Ulaby,  1981). 
The  effect  of  extinction  can  be  expressed  in  the  term  Ke  (extinction  coefficient).  Like 
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emissivity,  K:  is  a  value  between  0  and  1.  Because  it  depends  on  the  distance  traveled  in 
the  medium,  extinction  can  be  combined  with  distance,  z,  to  form  a  single  factor  which  is 
called  optical  depth,  t  where 

T(z,oo)  =  j^jc(z)^dz  (7) 

Note  from  Figure  6,  the  entire  atmosphere  affects  all  the  radiation  sources  from  the 
surface  (0)  to  the  top  of  the  atmosphere  (H).  The  optical  depth  is  sometimes  referred  to 
as  opacity  -  extent  to  which  a  given  layer  of  material  reduces  the  intensity  of  the  radiation 
(Rees,  1990).  The  atmosphere  reduces  in  density  and  hence  extinction  with  height.  The 
non-linear  decrease  of  extinction  can  be  expressed  with  an  exponential  function  as  e'^*’“\ 
This  reduction  factor  must  be  applied  to  all  sources  of  upwelling  radiation. 

The  surface  emission  (emissivity),  introduced  with  Equation  (5),  needs  expanding  to 
explain  medium  or  volume  emission.  Two  sources  of  radiation  affect  volume  emission. 
Note,  thermodynamic  equilibrium  is  assumed,  thereby  emission  equals  absorption.  The 
first  emission  source  is  the  radiation  emitted  by  the  matter.  The  second  source  to  enhance 
emission  is  scattering  of  radiation  into  the  direction  of  the  radiometer  from  other  matter 
outside  of  the  cylinder  viewed  by  the  radiometer  -  emission  scattering.  Like  opacity,  the 
upwelling  radiation  emitted  by  the  atmosphere  decreases  with  height.  In  contrast  to  the 
single  radiation  sources  at  the  surface,  a  sum  (integral)  of  sources  from  many  layers  is 
used  to  describe  the  atmospheric  emission.  The  atmosphere’s  emission  contribution  can 
be  expressed  using  the  product  of  K^and  its  brightness  for  each  increment  of  height,  dz. 
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However,  just  like  the  earth’s  emission,  the  atmosphere’s  emission  at  each  layer  is  also 
reduced  by  extinction  on  its  way  up  to  the  radiometer.  Thus,  the  reduction  factor, 
is  also  applied  to  the  atmospheric  source  term  as  it  was  to  the  surface  radiation  sources. 

In  essence,  the  reduction  factor  is  applied  twice.  However  for  microwave  radiometry,  the 
scattering  contribution  to  extinction  can  be  ignored.  As  Janssen  (1993)  states,  the 
common  remote  sensing  frequency  range  from  about  20  GHz  to  90  GHz  contains  certain 
frequencies  where  absorption  due  to  liquid  clouds  is  at  least  two  orders  of  magnitude 
greater  than  the  scattering.  If  scattering  is  neglected  error  rates  are  less  than  1%. 

By  using  Hollinger  (1983)  and  Ulaby  (1986)  as  a  guide,  the  above  effects  are  applied 
to  the  radiation  sources  in  a  radiative  transfer  equation.  Using  Equation  (4),  the  Rayleigh- 
Jeans  Law,  brightness  temperatures  represent  the  surface  and  atmosphere  source 
radiances.  The  equation  must  also  include  the  effects  of  emissivity  and  reflectivity.  The 
resulting  equation  is  rather  complex. 


H 

T,,  (H)  =  eTaC-^  +  +  J  Kj,_^^e-^dz  +  r 

0 

-0 

(8) 


where  Tap  =  Apparent  brightness  temperature  at  satellite  height  H  (K) 
e  =  Emissivity  of  ground 
Tg  =  Physical  temperature  of  ground  (K) 

e'^  =  Atmospheric  reduction  factor  (transmissivity  of  entire  atmosphere) 

T  =  Optical  depth  above  level  z 
r  =  Reflectivity  by  the  ground 

Tsky  =  Brightness  temperature  of  cosmic  and  galactic  sources 
Kc  =  Atmospheric  extinction  coefficient 

TA-down  &  Ta-up  =  Physical  temperature  of  infinitesimally  small  atmospheric  layer 
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Equation  (8)  represents  what  the  radiometer  sees  -  radiation  from  the  four  sources  in 
Figure  6.  The  first  term  is  the  surface-emission  blackbody  brightness  temperature  (equal 
to  the  physical  temperature)  reduced  by  its  emissivity  (non-blackbody  characteristic)  and 
the  intervening  atmosphere’s  transmission  (or  transmissivity).  The  second  term  is  the 
cosmic  and  galactic  background  radiation  reduced  three  times:  atmospheric  extinction 
inbound  to  the  earth,  reflection  off  the  earth,  and  then  atmospheric  extinction  again 
outbound  to  the  satellite.  The  third  term  represents  upward  atmospheric  emission.  The 
product,  K^TA-updz,  is  the  brightness  contribution  from  a  thin  atmospheric  layer  that  is 
reduced  on  the  way  up  to  the  radiometer  by  the  atmospheric  absorption  factor,  e'\  The 
integration  sums  the  brightness  contributions  for  all  atmospheric  layers.  The  last  term 
represents  downward  atmospheric  emissions  reduced  similarly  to  Tsky.  Note  all  the  optical 
path  terms  are  implicitly  multiplied  by  sec  0,  which  is  the  slant  path  up  to  the  radiometer. 
Equation  (8)  has  many  unknowns  and  the  intent  of  showing  the  reader  this  equation  is  to 
demonstrate  the  many  complexities  of  describing  the  radiation  viewed  by  a  satellite. 

The  radiative  transfer  equation  illustrated  by  Equation  (8)  shows  if  all  the  variables  on 
the  right-hand  side  of  the  equation  are  known,  then  the  brightness  temperature  measured 
at  the  satellite  can  be  calculated.  This  is  an  example  of  a  forward  problem.  In  reality,  the 
situation  is  an  inverse  problem.  The  satellite  provides  the  brightness  temperature  and  the 
unknowns  on  the  right  hand  side  of  Equation  (8)  must  be  determined  in  order  to  solve  for 
the  surface  temperature  (Janssen,  1993).  As  Ulaby  (1986)  points  out,  such  an  inverse 
problem  does  not  have  a  straightforward  solution. 
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The  treatment  of  the  unknowns  will  be  presented  in  a  later  section  describing  the 
SSM/I  surface  temperature  retrieval  technique.  Until  then,  exactly  what  the  satellite 
provides  as  output  must  be  explained.  Table  1  shows  some  typical  brightness 
temperatures  for  one  SSM/I  frequency  received  from  the  Air  Force  Weather  Agency 
(AFWA). 

In  reality  the  satellite  does  not  observe,  or  provide  as  an  output,  brightness 
temperatures.  The  next  section  briefly  outlines  the  principles  of  taking  radiometric 
measurements  and  how  the  brighmess  temperatures  in  Table  1  were  derived. 

Table  1  -  Brightness  Temperatures  for  the  19  GHz  SSM/I  Frequency 


Brightness  Temp.  (K) 

I-Coord. 

J-Coord. 

K- 

Coord. 

134.0 

1 

1 

43 

135.0 

2 

1 

43 

136.0 

3 

1 

43 

Antenna  Theory 

The  function  of  a  radiometer  is  to  measure  brightness  (or  radiant  intensity).  The 
brightness  is  an  idealized  flow  of  power  in  infinitesimally  small  elements  of  bandwidth 
(Janssen,  1993).  The  passive  antenna  can  be  defined  as  a  region  of  transition  between  a 
propagating  electromagnetic  wave  in  space  and  a  guided  wave  propagating  in  a 
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transmission  line  (voltage  in  a  wire)  (Ulaby,  1981).  Antennas  are  usually  made  of  metal  so 
as  to  conduct  electron  voltages.  The  electromagnetic  wave  induces  a  fluctuating  current 
in  the  antenna.  This  current  or  power,  P,  is  represented  by 

P  =  r„A,F.Wv  (9) 

where  Tap  =  Apparent  brighmess  temperature 
Af  =  Antenna  physical  receiving  area 
Fn  =  Antenna  electromagnetic  pattern 
k  =  Boltzmann’s  constant 
Av  =  Bandwidth 

As  in  Equation  (8),  Tap  is  the  apparent  brightness  temperature  measured  by  the  satellite 
radiometer.  Af,  Fn,  and  Av  are  known  from  the  engineering  of  the  antenna.  The  power  or 
voltage  counts  for  each  of  the  SSM/I  frequencies  is  the  actual  output  of  the  radiometer. 
These  voltage  counts  are  transmitted  to  various  sites  on  the  earth.  Using  the  voltage 
counts,  and  solving  for  Tap  in  an  equation  like  Equation  (9),  the  P  values  are  converted  to 
brighmess  temperatures  (as  in  Table  1)  at  AFWA,  Offutt  AFB. 

The  brightness  temperamres  are  then  mapped  onto  a  polar  stereographic  grid  (Figure 
7)  which  is  512  X  512  grid  points  in  dimension  for  a  hemisphere.  Each  box  has  a  grid 
dimension  of  64  x  64.  Note  there  are  two  grid-referencing  systems  for  Figure  7.  First,  a 
point  on  the  globe  can  be  represented  by  the  I,  J,  K  coordinates.  For  example,  a  point  in 
the  western  CONUS  might  be  1, 39, 44.  The  44  coordinate  corresponds  to  the  box  over 
the  central  US  and  Canada  as  shown  in  Figure  7.  A  point  could  also  be  represented  with 
just  two  coordinates,  I  and  J,  using  the  fact  the  whole  hemisphere  has  dimension  512  x 
512.  Thus,  the  same  point  mentioned  above  in  the  44  box  would  be  193, 359.  The 
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brightness  temperatures  mapped  onto  the  grid  are  referred  to  as  the  Sensor  Data  Record 
(SDR)  and  were  used  in  this  study  as  input  to  the  surface  temperature  algorithms. 

It  has  been  shown  what  the  radiometer  measures,  what  the  radiometer  provides  as 
output,  and  how  this  output  is  converted  to  a  brightness  temperature.  The  last  question  to 
answer  is  how  to  solve  the  inverse  problem,  i.e.,  how  to  obtain  surface  temperatures  from 
the  brightness  temperatures. 


Figure  7.  Polar  Stereographic  Northern  Hemisphere  Grid 
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Temperature  Retrieval  Technique  -  Statistical  Inversion  Method 

The  goal  is  to  find  a  function  that  relates  values  of  Tap  (v)  (the  measured  brightness 
temperatures  for  various  frequencies)  to  the  surface  temperature  for  various  land  types. 
The  function  can  be  written  in  a  manner  similar  to  Equation  (8)  where  the  brightness 
temperatures  are  written  on  the  right-hand  side  and  the  surface  temperatures  on  the  right. 
Writing  both  sets  of  temperatures  in  vector  form  yields: 


Tap  =  WT  (10) 

where  Tap  =  radiometric  observations  for  various  frequencies  from  a  particular  land  type 
W  =  weighting  matrix 

T  =  surface  temperatures  for  a  particular  land  type 

Tap  has  dimension  or  order  n  corresponding  to  the  number  of  microwave  frequencies, 

T  has  order  m  corresponding  to  the  number  of  surface  temperatures  taken  for  a  particular 
land  types,  and  W  is  of  dimension  mxn.  Equation  (10)  is  a  simple  representation  of  the 
complex  radiative  transfer  process  in  Equation  (8)  in  which  the  brighmess  temperatures 
derived  from  the  satellite  are  a  function  of  the  physical  temperature.  The  relation  or 
weighting  matrix  represents  all  the  other  radiation-matter  interaction  terms  in  Equation 
(8).  Equation  (10)  assumes  a  linear  relationship  and  Gaussian  distribution  of  temperatures 
(Ulaby,  1996).  The  Raleigh- Jeans  approximation  validates  the  linear  assumption.  The  fact 
that  atmospheric  parameters  follow  a  Gaussian  distribution  when  the  sample  size  is  large 
enough  to  evoke  the  central  limit  theorem  validates  the  Gaussian  assumption  (Wilks, 
1993).  Solving  for  the  desired  surface  temperatures  gives  Equation  (1 1). 
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T  =  W^Tap 


(11) 


There  are  many  inversion  techniques  that  provide  the  relation  matrix,  W*^  (Ulaby, 
1986).  The  SSM/I  program  uses  the  statistical  inversion  method  (SIM)  (HoUinger,  1983). 
SIM  essentially  determines  W‘  with  multiple-linear  regression.  Empirically,  several  sets 
of  radiometric  measurements  for  various  microwave  frequencies  (Tap)  is  obtained 
corresponding  to  simultaneous  measurements  of  T  (surface  temperatiures).  The  study 
obtains  several  sets  of  brightness  and  surface  temperatures,  each  set  corresponding  to  a 
different  land  type.  Standard  regression  programs  are  then  used  to  derive  the  relation 
matrix.  Hollinger  (1983)  used  a  very  similar  equation  to  Equation  (11)  for  application  to 
the  SSM/I, 


P  =  D  t  (12) 

where  P  =  vector  of  surface  temperature  estimates 
t  =  vector  of  brightness  temperatures 
D  =  correlation  coefficient  matrix  which  relates  P  &  t 

Note  there  is  no  actual  inversion  of  any  matrix  when  determining  the  brightness- 
physical  temperature  relationship.  “Inversion”  means  the  overall  problem  is  a  backward 
one  in  terms  of  standard  radiative-transfer  equations  —  solving  for  the  right-hand  side 
versus  the  left-hand  side  of  the  equation.  Multiple-linear  regression  is  used  to  define  the 
relationship  between  brightness  temperatures  and  physical  temperatures.  The  regression 
coefficients  and  frequency  variables  are  then  put  in  algorithm  form  (as  shown  in  the  next 
chapter).  The  result  is  a  set  of  algorithms  where  each  algorithms  is  used  for  a  particular 
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land  type.  Operationally,  the  real  time  brightness  temperature  measurements  (t)  serve  as 
input  to  the  algorithms. 

Janssen  (1993)  adds  an  error  term  (e  =  £1+62+63...)  to  equation  (12).  The  errors  could 
be  caused  by  several  reasons  inherent  to  remote  sensing:  radiometer  noise,  linearization, 
unknown  parameters,  and  weighting  functions.  The  end  of  this  chapter  contains  a  fuU 
discussion  of  errors. 

With  the  generation  of  the  algorithm  set  complete,  where  each  algorithm  pertains  to  a 
particular  land  type,  the  land  type  must  be  known  to  retrieve  a  surface  temperature.  An 
empirical  approach  is  used  to  determine  the  land  type  relationship.  The  power  or 
brightness  temperature  from  each  channel  usually  falls  within  a  certain  range.  These 
ranges  are  experimentally  determined  for  various  land  types  and  then  combinations  of 
these  ranges  are  put  in  mle  form.  For  example,  the  land  type  dry  arable  soil  is  detenmned 
in  the  CV  algorithm  with  the  following  rule: 


If  (V22  -  V19)  <  4.0  and 
4.0  <  (V19  +  V37)/2  <  9.8  and 
(V37  -  V19)  >  -6.5  and 
-5.0  <  (V85  -  V37)  <  0.5  and 
(V85  -  H37)  <  4.2  then 
land  type  =  dry  arable  soil 

where  V22  =  Vertically  polarized  22.2  GHz 
V19  =  Vertically  polarized  19.3  GHz 
V85  =  Vertically  polarized  85.5  GHz 
V37  =  Vertically  polarized  37.0  GHz 
H37  =  Horizontally  polarized  37.0  GHz 

With  similar  rules  established  for  the  various  land  types,  the  rules  are  applied  to  the 
SSM/I  channel  measurements  to  determine  the  land  type.  With  the  land  type  determined. 
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the  corresponding  surface  temperature  algorithm  is  used  to  calculate  the  surface 
temperature.  The  land  type  rules  are  not  all  inclusive.  Some  channel  measurements  do 
not  fall  into  one  of  the  land  type  rales  and  thus  a  land  type,  and  in  turn,  a  surface 
temperature  is  undeterminable.  As  a  result,  the  land  typing  schemes  affect  the  quantity  of 
algorithm-produced  surface  temperatures. 

In  summary,  the  method  used  with  SSM/I  data  to  calculate  surface  temperatures  does 
not  require  prior  (a  priori)  knowledge  of  the  unknowns  in  Equation  (8):  temperature 
profiles,  emissivity,  reflectivity,  absorption,  scattering,  etc.  of  the  earth-atmosphere 
system.  Instead,  what  is  needed  is  a  priori  knowledge  of  the  relation  between  measured 
brighmess  temperature  at  seven  channels  and  measured  physical  temperature  for  various 
land  types.  Multiple  linear  regression  determines  those  relationships.  The  regression 
coefficients  are  used  to  build  sets  of  TS  and  CV  surface  temperature  algorithms,  each 
algorithm  corresponding  to  a  particular  land  type.  A  priori  knowledge  of  the  land  typing 
schemes  is  required  to  choose  which  algorithm  to  use  in  the  set. 

Calibration/Validation  (CV)  Algorithm  and  FORTRAN  Code 

Still  in  wide  use,  the  CV  algorithm  in  this  study  is  based  upon  the  land  surface 
temperature  algorithm  as  outlined  by  McFarland  (1991).  The  Naval  Research  Laboratory 
(NRL)  produced  CV  in  1988  to  calibrate  and  validate  the  SSM/I.  The  algorithm  is  a 
linear  combination  of  coefficients  applied  to  the  SSM/I  frequencies,  but  only  4  channels 
are  used:  19H,  22V,  37V,  and  85V.  Different  algorithms  are  used  based  upon  one  of  four 
different  land  types  classifications  as  shown  in  Table  2. 
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The  CV  program  was  quite  concise  -  only  four  pages  of  code.  The  program  consisted 
of  two  main  routines.  First,  the  land  type  was  determined  using  eight  land  type  rules,  as 
shown  previously  for  dry  arable  soil.  The  similar  land  type  classifications  (e.g.  moist  soils) 
were  combined  into  the  four  general  land  type  classifications  listed  in  Table  2.  Second, 
based  upon  one  of  the  four  general  land  types,  the  corresponding  algorithm  from  Table  2 
was  run  to  obtain  the  surface  temperature.  As  mentioned,  the  land  type  classifications 
and/or  rules  were  not  all  inclusive.  Occasionally,  the  program  could  not  determine  a  land 
type,  and  in  turn  a  surface  temperature  could  not  be  calculated. 

Moncet  et  al  (1996)  states  the  major  advantages  of  the  relatively  simple  statistical 
inversion  method  for  the  SSM/I  derived  parameters  are  robustness  and  operational 
expediency.  The  CV  code  is  certainly  expedient.  However,  he  points  out  the  dependence 
on  a  linear  relationship  and  lack  of  a  non-linear  treatment  of  the  background  conditions 

Table  2.  CV  Land  Surface  Temperature  Algorithms 


LAND  TYPE 

ALGORITHM 

Dense  Vegetation 

T(K)  =  24.94  -  1.2784(T19H)  +  .88(T22V)  +  .5933(T37V)  +  .7299(T85V) 

Agricultural/Range 

T(K)=  6.97  -  .6266(T19H)  +  .2716(T22V)  -  .1297(T37V)  +  1.482(T85V) 

Moist  Soils 

T(K)  =  23.16  -  .1873(T19H)  +  .5221(T22V)  -  .6271(T37V)  + 

1.232(T85V) 

Dry  Soils 

T(K)  =  72.68  -  .4598(T19H)  +  .5984(T22V)  +  .8828(T37V)  - 

.2623(T85V) 
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(e.g.  emissivity,  reflectivity)  may  lead  to  a  large  spread  of  retrieval  errors  (see  error 
discussion  at  end  of  chapter).  McFarland  (1991)  points  out  the  downside  of  this  technique 
for  retrieving  land  surface  temperatures  as  the  simple  treatment  of  the  land  types.  Land 
surfaces  not  only  vary  a  great  deal  spatiaUy,  but  also  temporally.  Temporal  variances 
occur  throughout  the  year  and  from  year  to  year  both  from  natural  and  man-made  changes 
in  the  land  surface. 

Temperature  from  Sensor  Microwave  Imager  (TS)  Algorithm  and  FORTRAN  Code 
In  contrast  to  the  CV  algorithm,  the  TS  algorithm  uses  all  seven  channels  to  obtain  a 
surface  temperature  (See  Table  3  below).  Additionally,  TS  appUes  eight  land  types  with 
the  potential  for  a  total  of  25  different  land  surfaces.  Currently,  the  snow  and  ice 
background  land  types  are  turned  off  due  to  unreliability  (Kopp,  1997).  Note  the  TS 
algorithms  incorporate  the  viewing  angle  (0)  as  well.  Another  difference  between  CV  and 
TS  is  AFWA  calculates  new  TS  coefficients  for  each  new  DMSP  satellite.  Theoretically 
each  new  satellite  and  its  SSM/I  has  a  slightly  different  viewing  angle  from  the  previous 
satelUte  and  thus,  the  relation  between  brightness  temperature  and  physical  temperature 
should  change  as  well.  CV  coefficients  and  their  algorithms  have  been  constant  since  their 
inception  in  1987. 

Three  times  the  length  of  CV,  the  TS  program  (15  pages)  is  significantly  more  complex 
than  the  CV  program.  However,  the  TS  code  is  tunable  which  allows  for  continuous 
improvement  over  time.  CV  is  not  tunable. 
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Table  3.  TS  Land  Surface  Temperature  Algorithms  for  DMSP  F13  (provided  by  AFWA) 


LAND  TYPE 

ALGORITHM 

Densely  &  Moderately 

Vegetated  Land,  Cold 

Vegetated  Land 

T(K)=  18.833+  .8557(T22V)  -  .23(T37V)  +  .9349(T85V) - 

.6188(T85H)  +.4235(0) 

Arable  Soil 

T(K)  =  29.926  -  1.087 (T19H)  +  3.276(T22V)  -  1.664(T37H)  + 

10.357(T85H)  -  1.829(T85V)  +  5.243(0) 

Snow  &  Deep  Snow 

Covered  Land 

Not  Used 

Glacial  Snow 

Not  Used 

Morning  Desert 

T(K)  =  133.173  +  2.07(T19H)  +  7.464(T22V)  -  9.741(T37V) 

+  2.016(T37H)+6.967(T85V)  -  3.155(T85H)  +5.324(0) 

Sparsely  Vegetated 

Land 

T(K)  =  73.304  +  7.456(T22V)  -  7.546(T37V)  +  2.557(T35H) 

+9.275(T85V)  -  3.993(T85H)  +1.431(0) 

Flooded  Land 

T(K)  =  -4.301  -  2.554(T19H)  +  12.837(T22V)  -  8.893(T37V) 

+  2.889(T37H)  +  14.717(T85V)  -  8.805(T85H)  +4.555(0) 

Desert,  Wet  Surface 

T(K)  =  109.079  -  7.194(T19V)  +  .2841(T19H)  +  .6520(T22V) 

-  2.919(T37V)  +  14.223(T85V)  -  .6845(T85H) 
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Applications  of  SSM/I  Data 

This  paper  addresses  two  applications  of  the  algorithms  in  question,  both  of  which  are 
currently  operational  at  the  AF  Weather  Agency  (AFWA).  The  TS  algorithm  is  primarily 
used  in  two  capacities.  First,  TS  serves  as  one  input  for  the  global  surface  temperature 
model.  Surface  Temperature  (SFCTMP).  Secondly,  TS  aids  the  analysis  of  global  cloud 
coverage  (Kopp,  1995).  Existing  as  a  sub-set  of  a  cloud  model  called  Real-Time 
Nephanalysis  (RTNEPH),  TS  was  rewritten  in  1990  as  part  of  an  overall  improvement 
effort  to  analyze  and  forecast  clouds  (Kopp,  1995).  To  analyze  clouds  in  data  sparse 
areas,  AFWA  uses  the  TS  output  to  first  determine  if  clouds  are  present  by  comparing  the 
IR  temperature  to  the  background  surface  temperature  partly  derived  from  SSM/I  data.  If 
the  IR  temperature  is  colder  than  the  surface  temperature  then  a  cloud  is  assumed  to  be 
present  (Kopp,  1994).  If  the  derived  background  temperature  is  the  same  as  the  IR 
temperature,  when  in  fact  the  real  background  surface  temperature  is  warmer  than  the  IR 
temperature,  the  analysis  will  falsely  show  no  cloud  present. 

The  second  application  addressed  here  is  the  desire  of  military  forecasters  to  have  near 
real-time  surface  temperatures  from  data  sparse  or  denied  territory.  Among  many 
applications,  surface  temperature  is  a  key  element  in  the  ability  of  precision  guided 
munitions  to  acquire  their  targets.  Heavily  based  upon  the  temperature,  forecasters  wiU 
brief  if  there  is  enough  thermal  contrast  between  the  target  and  it’s  background.  To  assist 
in  the  determination  of  surface  temperatures,  the  USAF  is  currently  fielding  its  newest 
deployable  satellite  imagery  receiver  called  the  Small  Tactical  Terminal  (STT).  The  STT 
uses  the  CV  algorithm  to  create  the  surface  temperature  products  (Harris  Corp.,  1997). 
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Deployed  forecasters  need  to  know  the  average  departures  from  truth  (bias)  and  accuracy 
of  the  surface  temperatures  produced  by  CV. 

Sources  of  Error 

According  to  Janssen  (1993),  radiometric  errors  can  arise  from  essentially  four 
sources:  radiometer  noise,  linearization,  unknown  parameters,  and  regression. 

The  first  is  radiometer  noise,  which  is  composed  of  random  and  bias  components. 
Random  noise  is  Gaussian  in  nature  and  thus  has  known  statistics.  However,  bias  is  not 
obvious  and  must  be  reduced  to  tolerable  levels  by  calibration  and  comparison  of  the 
derived  products  with  ground  truth  data.  Bias  will  be  uncovered  in  this  study.  AFWA 
regularly  calculates  and  removes  the  bias. 

Second,  linearization  errors  may  occur  due  to  the  Raleigh- Jeans  approximation.  If 
these  errors  are  not  less  than  the  radiometer  noise  then  other  inversion  methods  must  be 
applied. 

Third,  errors  from  unknown  parameters  exist.  For  example,  the  land  typing  is  not 
exact  due  to  the  high  spatial  and  temporal  variability  of  land  types.  Another  example  of 
errors  from  unknown  parameters  is  changes  in  satellite  viewing  angle.  Fuhrhop  and 
Simmer  (1995)  performed  a  study  on  SSM/I  radiometers  and  found  the  variations  of 
viewing  angle  to  be  as  high  as  1.5  degrees  which  results  in  a  brightness  temperature 
change  of  2  K.  By  analyzing  the  brightness  temperature  variations  with  changing 
incidence  angle,  a  correction  can  be  derived  and  used.  TS  corrects  for  variable  viewing 
angle  by  including  it  in  the  regression.  A  third  example  of  unknown  parameters  error  is 
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the  effect  of  hydrometers  -  specifically  cloud  water  droplets  and  rain.  As  Ulaby  (1981) 
points  out,  water  clouds  reduce  microwave  transmission  below  wavelengths  of  2  cm,  thus 
affecting  all  the  channels  on  the  SSM/I  (wavelengths  given  below  Figure  6).  Ice  clouds 
have  almost  no  effect.  A  major  source  of  error  is  rain.  The  rain  attenuation  is  appreciable 
at  wavelengths  below  2  cm.  However  the  amount  of  intense  rain  actually  covering  the 
earth  to  affect  SSMI  readings  is  quite  small  (Ulaby,  1981).  Janssen  (1993)  points  out 
taking  measurements  at  other  frequencies  reduce  the  errors  from  unknown  parameters. 

Finally,  there  are  potential  errors  in  calculation  of  the  regression  coefficients.  The 
multiple  linear  regression  is  only  an  approximation  to  the  brightness  temperature-physical 
temperature  relationship.  The  regression  analysis  has  some  associated  error. 

Remote  Sensing  of  Land  Surface  Temperatures  -  Literature  Review 

The  science  of  Passive  Microwave  Radiometry  (PMR)  has  not  focused  on  the  retrieval 
of  surface  temperatures  over  land  in  contrast  to  ocean  parameters  and  atmospheric 
profiles.  Indeed,  there  is  no  hard  DoD  requirement  for  the  performance  of  the  SSM/I 
derived  surface  temperatures  (Glackin  and  Gasster,  1994).  However,  within  the  last 
several  years  there  have  been  some  related  studies. 

Over  the  continent  of  Australia,  Prata  (1994)  investigated  PMR  retrieval  of  land 
surface  temperatures,  but  his  study  used  data  from  NOAA’s  Advanced  Very  High 
Resolution  Radiometer  (AVHRR).  The  AVHRR  uses  the  IR  portion  of  the  EM  spectrum, 
not  the  microwave  portion.  Also  using  IR,  but  from  radiometers  onboard  aircraft,  Chen 
and  Engman  (1997)  studied  the  spatial  pattern  of  moisture  and  temperature  over  the 
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central  US  prairie.  They  show  their  algorithms  to  generally  overestimate  the  surface 
temperature  by  2.5  degrees  with  an  RMSE  of  3.51  K.  A  study  by  Macelloni  et  al  (1994) 
used  SSM/I  data  to  determine  land  types  but  not  surface  temperatures  in  several  test  areas 
around  the  globe.  The  original  SSM/I  calibration  and  validation  (1987)  of  the  land  surface 
temperature  algorithm  is  the  study  most  similar  to  this  thesis.  Although  using  DMSP 
satellite  F8  and  covering  only  over  the  Western  Desert  and  Central  Plains  of  the  US,  the 
Naval  Research  Lab  collected  a  total  of  1764  data  pairs  for  4  days  in  Aug  1987.  The  CV 
algorithm  performed  with  an  RMSE  of  3.1  K  for  all  land  types  (McFarland,  1991).  Thus, 
the  initial  work  by  Comoglio  (1997)  and  continued  work  by  this  author  appear  to  be 
unique. 

Summary 

For  a  given  frequency,  radiative  transfer  theory  shows  the  physical  temperature  is 
related  to  brightness,  and  that  relationship  is  linear  for  the  microwave  portion  of  the  EM 
spectrum  at  Earth  temperatures.  For  sensing  land  surface  temperatures  from  space,  the 
brightness-physical  temperature  relationship  is  complex  when  specifying  all  the  sources 
and  radiation-matter  effects. 

The  SSM/I  project  overcomes  this  complex  relationship  empirically.  Multiple  linear 
regression  applied  to  field  measurements  determines  the  temperature-brightness 
relationship  in  a  set  of  algorithms  for  various  land  types.  Operational  microwave 
brightnesses  serve  as  input  to  the  algorithms,  which  then  provide  the  surface  temperatures. 
The  question  remains  as  to  the  performance  level  of  the  TS  and  CV  algorithms. 
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in.  Methodology 


Chapter  Overview 

The  methodology  section  explains  the  processing  of  the  data  and  the  statistics  gathered 
from  the  data.  The  main  goal  of  data  processing  was  to  obtain  data  pairs  -  ground  station 
and  algorithm  temperature  at  the  same  time  and  place.  Not  withstanding  inherent  errors, 
ground  station  data  was  assumed  to  be  "truth"  data. 

The  following  chapter  starts  with  how  the  current  work  fits  in  with  prior  AFIT 
research.  Data  sources  are  then  described  followed  by  a  detailed  explanation  of  the  data 
comparison  process  that  generates  data  pairs.  The  chapter  ends  with  a  discussion  of 
algorithm  performance  statistics. 

Building  on  Past  Work  at  AFIT 

Comoglio’s  (1997)  work  was  essentially  a  preliminary  study  of  algorithm  performance 
because  most  of  his  work  entailed  the  proper  setup  and  execution  of  the  algorithms  on 
AFTT  workstations.  Time  was  not  available  to  accomplish  an  extensive  comparison  over 
time  and  space.  The  task  of  porting  the  algorithms  to  AFIT  hardware  alone  was 
considerable,  especially  for  the  TS  algorithm.  Since  the  present  work  was  a  continuation 
of  Comoglio’s  research,  it  was  imperative  to  precisely  follow  his  procedures  to  properly 
combine  the  results. 
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To  prove  continuity  of  Comoglio’s  work,  the  CONUS  summer  96  data  set  was  chosen 
to  be  a  trial  run  of  the  process.  Table  4  displays  both  the  Comoglio  and  current  study’s 
CONUS  Aug  96  results.  The  process  was  duplicated  well. 

This  study  accomplished  major  improvements  in  processing  and  data  organization. 
Instead  of  manually  entering  each  hour  of  a  day,  programs  were  made  to  automatically 
cycle  through  all  24  hours.  In  addition,  automation  was  applied  to  the  SSM/I  data  files. 
Before  each  file  pertaining  to  a  particular  frequency  was  processed  separately.  During  the 
current  study,  all  SSM/I  files  for  one  day  were  automatically  processed  together. 

Data  Sources  and  Description 

The  SSM/I  data  for  the  Northern  Hemisphere  was  recorded  on  8mm  tapes  at  AFWA, 
Models  Branch,  Offutt  AFB  Nebraska  and  sent  to  AFIT.  AFWA  provided  the  data  in 
UNIX  gridded  format  as  described  in  Chapter  2. 

The  SSM/I  data  set  was  organized  by  day.  Each  day  had  3  file  types  comprising  a  total 
of  9  files  per  day:  grid  date  time  stamp;  brighmess  temperatures  -  7  files,  one  for  each 
channel;  and  satellite  identifier.  See  Table  5. 
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Table  5.  SSM/I  Raw  Data  Files 


FILE  TYPE 

FILE  NAME 

Grid  Date  Time  Stamp 

RNXMI1_00MITT 

19  GHz  -  Vertical 

RNXMIl.OOMIHl 

19  GHz  -  Horizontal 

RNXMI1_00MIV1 

22  GHz  -  Vertical 

RNXMIl_00MrV2 

37  GHz  -  Horizontal 

RNXMI1_00MIH3 

37  GHz  -  Vertical 

RNXMI1_00MIV3 

85  GHz  -  Horizontal 

RNXMI1_00MIH4 

85  GHz  -  Vertical 

RNXMI1_00MIV4 

Satellite  Identifier 

RNXMI1_00MIID 

The  AF  Combat  Climatology  Center  (AFCCC),  Asheville,  North  Carolina  provided  the 
ground  station  data  in  DOS  format  as  shown  in  Table  6. 


Table  6.  Ground  Station  Data  Formats 


1— block# — 1 1 — Date - 1  l-time-l  l-tempK-l 

Raw 

6901901996080100003042 

Decoded 

690190, 1  Aug  96, 00  hr,  00  min,  304.2  K 

The  block  #  refers  to  the  World  Meteorological  Organization’s  geographic  reference 
system.  AFCCC  also  provided  soft-copy  of  station  identification  tables  which  were  used 
to  convert  the  station  block  number  to  the  more  familiar  three  letter  station  identifier  or 
International  Civil  Aviation  Organization  (ICAO)  code  name.  Data  was  sent  on  3.5-inch 
disks  in  compressed  format. 
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The  Data  Comparison  Process 

The  data  comparison  process  involved  manipulating  SSM/I  and  ground  station  or 
“truth”  databases  in  order  to  obtain  paired  data.  Again,  a  data  pair  consists  of  truth  and 
algorithm  surface  temperatures  at  the  same  time  and  place.  The  process  can  be  visualized 
in  Figure  8.  The  process  consisted  of  five  phases:  Download,  Reformat,  Calculate,  Match, 
and  Compare.  The  first  three  phases  were  applied  to  both  the  raw  SSM/I  data  and  the 
ground  station  data.  The  fourth  phase  created  the  data  pairs  with  a  matching  process. 

The  fifth  phase.  Compare,  calculates  the  statistics  from  the  data  pairs. 

In  terms  of  software,  phases  two,  three,  and  four  required  a  FORTRAN  program  to 
execute  (see  Figure  8).  Eight  FORTRAN  programs  were  created  to  complete  the  data 
processing.  With  the  exception  of  the  algorithms,  Comoglio  and  this  author  created  all  of 
the  code  for  the  data  comparison  process.  Sun  Sparc  workstations  provided  the  platform 
for  phases  1-4.  The  comparison  phase  occurred  on  a  PC  using  standard  software. 

Detailing  each  phase,  step  A  consisted  of  transferring  the  raw  SSM/I  data  on  8mm  tape 
and  the  station  data  on  3.5  diskette  to  the  workstations.  The  SSM/I  data  was  transferred 
from  tape  directly  to  the  workstations  and  a  working  directory  using  the  UNIX  utility 
program  tar  (Sobel,  1985).  The  station  data  was  first  uncompressed  from  a  PC  on  the 
AFIT  LAN,  then  transferred  to  the  UNIX  workstations  using  file  transfer  protocol. 
Station  files  were  converted  to  UNIX  format.  The  next  phase  in  the  process  was 
reformatting. 

The  raw  SSM/I  data  needed  reformatting  from  binary  to  ASCII.  In  step  B,  a  program 
(SATTIME.F)  filtered  the  satellite  date- time  data  for  the  area  of  interest  (grid  box)  and 
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4.  Match 

5.  Compare 


5  Phases 

1.  Download 

2.  Reformat 

3.  Calculate 


Match  TB  by  I,  J,  K  w/ MATCH.F 

Stn  T  TS  T 

CV  T 

IS-Neph 

J  S-Nqph 

DYS  280.4  273.2 

274.5 

211.04 

386.73 

CHS  284.3  277.3 

278.9 

236.4 

375.5 

Figure  8  -  Data  Comparison  Process 


day.  The  output  is  in  ASCII  format  in  a  separate  file  in  the  form  shown  in  Table  7. 


Table  7.  Output  from  Satellite  Observation  Time  Program 


JuUan 

Day  YRHRMin 

I  J  K 
coordinates 

214960307 

1  1  43 

214960307 

2  1  43  ... 

In  step  C,  another  program  (SSMI.F)  filtered  the  channel  data  files  for  area  and  day  of 
interest,  added  the  decimal  point  to  the  brightness  temperature  (see  1**  column  in  Table  8), 
and  reformatted  the  data  to  an  ASCII  file.  Steps  A,  B,  and  C  were  complete  up  to  this 
point.  Steps  D  and  E,  the  surface  temperature  calculations,  were  next. 

Phase  3,  the  calculation  part  of  the  process,  consisted  of  the  actual  algorithms  that 
produce  the  surface  temperatures.  The  algorithm  programs  produced  output  consisting  of 
temperatures  for  each  grid  point  such  as  shown  in  Table  8. 


Table  8.  Output  from  Surface  Temperature  Algorithm  (CV  or  TS)  Programs 


Temp  (K) 

I  J  K 
coordinates 

134.0 

1  1  43 

135.0 

2  1  43... 

With  phase  3  complete,  programs  were  run  to  calculate  the  production  rates  for  each 
algorithm.  The  production  rate  represented  the  percent  of  all  the  grid  points  in  the  area  of 
interest  which  returned  a  realistic  temperature  (step  F).  By  this  point,  all  the  SSM/I  data 
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had  been  processed  and  was  ready  for  matching  (Phase  4).  Next,  the  station  data  needed 
processing  to  prepare  for  the  matching  phase. 

After  downloading  and  transforming  from  diskettes,  the  station  data  was  filtered  by 
hour  of  interest  corresponding  to  the  satellite  observation  times.  Step  H  was 
accomplished  by  the  program  called  BLOCKNUM.F.  For  one  station  observation  at  a 
time,  the  program  converted  the  temperature  and  latitudeAongitude  to  real  numbers  by 
adding  the  decimal  point,  converted  latitude/longitude  to  a  polar  stereographic  grid 
coordinate,  and  created  a  new  data  file  ready  to  compare  to  algorithm  temperature  files. 
Output  had  the  following  items:  ICAO,  temperature  (K),  latitude,  longitude,  I-J-K 
coordinates  given  in  the  form  as  shown  in  Table  9. 


Table  9.  Ouq)ut  from  Ground  Station  Program 


ICAO 

Temp  (K) 

Latitude  (°) 

Longitude  (°) 

I-Coord. 

J-  Coord. 

DYS 

299.2 

32.26 

99.51 

211.04 

386.73 

CHS 

299.9 

32.54 

80.02  1 

256.95 

393.84... 

Phase  4  entailed  the  critical  matching  phase.  Creating  the  algorithm-station 
temperature  data  pairs  was  done  in  three  general  routines  by  the  matching  program 


(MATCH.F): 

1  _  Por  a  particular  day  and  hour,  SSM/I  temperatures  and  location  were  obtained  by 
matching  the  coordinates  of  the  satellite  observation  time  with  the  algorithm- 
produced  temperature  and  grid  coordinate  (files  from  steps  B  &  C  are  combined). 
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2.  The  station  observations  are  filtered  for  hour(s)  of  interest  to  obtain  the 
corresponding  temperature  and  grid  coordinate. 

3.  Lastly,  the  algorithm  surface  temperature  and  station  surface  temperature  are 
matched  by  grid  coordinate. 


A  match  is  made  if  a  grid  with  a  temperature  is  within  one  grid  coordinate  of  a  ground 
station.  Reference  Figure  9  taken  from  Comoglio  (1997).  Thus,  for  each  of  the  ground 
stations  the  program  gathers  a  maximum  of  four  surrounding  points.  The  matching 
program  output  looks  like  Table  10. 


Table  10.  Output  from  Matching  Program 
showing  station,  TS,  and  CV  temperatures. 


ICAO 

Temp 

(K) 

TS 

Temp. 

(K) 

CV 

Temp. 

(K) 

I-Super 

Coord. 

J-Super 

Coord. 

I- 

Coord. 

J- 

Coord. 

DYS 

280.4 

273.2 

274.5 

211.04 

386.73 

248 

373 

CHS 

284.3 

277.3  ^ 

278.9 

236.43  ^ 

375.52~^ 

237 

376 

A  match  can  be  of  a  station  temperature  with  just  a  TS  temperature  or  a  CV 
temperature  or  both.  The  match  step  produced  hourly  files  by  day,  which  were  then 
manually  combined  into  one  file  using  a  text  editor.  Finally,  the  single  file  was  converted 
to  ASCII  format  for  import  into  various  PC  programs.  The  match  step  at  this  point  was 
completed. 
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After  the  combined  DOS  files  were  imported  into  Statistix  or  MS-Excel,  various 
descriptive  statistics  were  calculated.  At  this  point,  data  processing  was  complete.  The 
next  task  was  the  generation  of  the  statistics. 


Comparison  Methods 

This  section  outlines  what  was  done  to  the  large  set  of  data  pairs.  With  a  focus  on 
customer  and  operational  needs,  three  statistics  are  outlined:  production  rate,  bias,  and 
accuracy.  These  three  statistics  clearly  show  algorithm  differences.  Lastly,  data 
stratification  or  grouping  is  addressed. 
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Statistics 


The  appropriate  statistics  were  mostly  chosen  based  upon  customer,  i.e.,  AFWA, 
requirements.  As  mentioned,  the  DoD  does  not  have  formal  performance  requirements 
for  SSM/I  land  surface  temperatures  (Glackin  andGasster,  1994). 

The  first  descriptive  statistic  generated  was  production.  This  answered  the  first  major 
question  of  the  research,  “How  much  data  do  the  algorithms  produce?”  The  algorithms 
did  not  produce  temperatures  at  every  grid  point  due  to  the  inability  to  detenmne  the  land 
type.  Production  was  simply  the  percent  of  grid  points  with  a  temperature. 

The  second  statistic  generated  was  bias.  The  weather  forecaster  using  a  SSM/I  surface 
temperature  product  needs  to  know  the  tendency  or  bias  of  the  algorithms.  That  is,  are 
the  retrieved  temperatures  usually  too  cold,  too  warm,  or  equal  to  the  actual  surface 
temperature?  AFWA  did  not  require  the  bias  statistics  because  bias  is  calculated  and 
corrected  for  at  AFWA.  Also  called  mean  errors  (Wilks,  1995),  biases  were  calculated 
using  the  following  formula  fi*om  Kopp  (1994): 

Bias  =  ’^^AT  (15) 

N 

where  AT  =  Algorithm  temperature  minus  station  temperature 
N  =  Number  of  data  pairs  for  the  particular  algorithm 

The  last  statistic  generated  was  accuracy.  Most  importantly  for  AFWA,  accuracy 
answers  the  second  question,  “How  good  are  the  algorithms?”  As  opposed  to  production, 
accuracy  meant  quality  of  the  data.  The  Root  Mean-Square  Error  (RMSE)  defines 
accuracy  (Kopp  1994). 
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Accuracy  = 


(16) 


Accuracy  referred  to  the  average  correspondence  between  individual  algorithm 
temperature  observations  and  the  events  they  simulate  (truth  or  station  data).  In  contrast, 
bias  (or  systematic  bias)  measured  the  correspondence  between  average  algorithm 
temperature  and  the  average  station  temperature.  The  surface  temperature  algorithms 
which  were  consistently  too  warm  or  cold  exhibit  a  bias,  regardless  of  the  algorithm’s 
accuracy  (Wilks,  1995).  Of  important  note,  AFW A  stated  an  RMSE  >  8  K  was 
unacceptable.  Both  algorithms  were  compared  to  AFWA’s  accuracy  criteria. 


Data  Stratification 

Since  the  data  covered  a  time  span  of  at  least  seven  months  and  had  global  distribution, 
the  data  could  be  sorted  or  stratified  either  in  time  or  space.  Recall  SSM/I  data  was 
available  for  the  entire  Northern  Hemisphere.  Spatially,  the  data  was  grouped  into  four 
regions.  The  first  data  set  was  over  the  Continental  United  States  (CONUS).  The  study 
was  then  expanded  to  include  specific  areas  with  particular  interest  to  the  United  States 
Department  of  Defense:  Korean  Peninsula,  Saudi  Arabia,  and  Bosnia.  Descriptive  statistics 
were  calculated  for  three  seasons;  summer,  fall,  and  winter  for  each  region. 

Unfortunately,  the  temperature  data  was  too  sparse  to  provide  a  color  enhanced 
pictorial  representation  of  brightness  temperatures  for  the  satellite  passes.  In  addition,  not 
all  matches  for  one  location  contained  pairs  with  both  CV  &  TS  temperatures.  The 
matches  themselves  were  very  sparse  in  time  (see  appendices).  Thus,  comparative 
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statistics  could  not  be  shown  pictorially  as  well.  A  temporal  stratification  was  envisioned 
such  as  breaking  the  seasons  up  into  weekly  or  daily  sets.  However,  time  did  not  perimt 
such  stratification,  nor  was  AFWA  concerned  with  one.  The  primarily  concern  was 
seasonal  and  regional  algorithm  performance. 

Summary 

A  complex  process  with  many  separate  filter  and  reformatting  programs  was  required 
to  obtain  the  data  pairs  of  truth  and  algorithm  surface  temperature.  Once  the  data  pairs 
were  gathered,  the  customer  driven  statistics  were  calculated  such  ^  production,  bias,  and 
RMSE  values.  With  all  the  data  processed  and  statistics  gathered  the  analysis  process 
could  begin. 
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IV.  Results  Analysis 


Chapter  Overview 

This  chapter  has  five  main  sections.  After  a  short  discussion  about  the  data  sets,  the 
algorithm  production  is  reviewed.  The  next  section  discusses  the  tally  of  data  pairs.  The 
last  two  sections  reveal  the  variability  of  the  algorithm  bias  and  accuracy. 

The  sections  with  production,  bias,  and  accuracy  results  follow  a  basic  pattern.  The 
first  sub-section  provides  the  overall  results,  e.g.,  grand  averages.  Then,  the  algorithm 
numerical  differences  as  compared  to  each  other  are  detailed.  The  last  sub-section 
discusses  any  regional  and/or  seasonal  patterns. 

Data  Sets 

Twelve  data  sets  were  created  -  four  regions  each  with  three  seasons.  The  four 
regions  were  the  CONUS,  Bosnia,  Korean  Peninsula,  and  Saudi  Arabia.  In  reality,  Bosnia 
was  composed  of  some  surrounding  countries,  which  were  part  of  the  former  Yugoslavia. 

The  number  of  data  days  for  each  season  varied  from  region  to  region  depending  upon 
SSM/I  data  availability.  In  general,  summer  (1996)  had  14  days,  fall  (1996)  had  15  days, 
and  winter  (1996-1997)  had  18  days  (reference  appendices). 

Algorithm  Production 

The  first  research  question  addressed  was,  which  algorithm  produced  more 
temperatures?  Production  is  defined  as  the  percent  of  temperatures  calculated  in  the  grid 
box  (s)  of  interest.  The  important  factor  is  the  relative,  i.e.,  TS  vs  CV,  production  rate. 
Absolute  production  rate  (out  of  100)  is  irrelevant  as  the  grid  boxes  usually  overlay  some 
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ocean  areas,  which  naturally  reduce  the  number  of  land  surface  temperatures.  The 
production  rates  are  found  in  Table  11.  Production  highlights  only  algorithm  temperature 
generation  and  not  any  comparison  with  ground  observations. 


Table  11.  Surface  Temperature  Production  Rate  (%)  for  the  Temperature  from  Sensor 
Microwave  Imager  (TS)  and  Calibration/V alidation  (CV)  Algorithms. 

Each  region  shows  its  approximate  central  latitude 


SUM 

MER 

FALL  1 

WINTER 

Reqion  &  Latitude 

TS 

CV 

TS 

CV 

TS 

CV 

33 

17 

28 

28 

18 

21 

70 

51 

69 

55 

21 

32 

KOREA  (37°) 

84 

50 

40 

41 

8 

26 

SAUDI  ARABIA 
(24°) 

81 

63 

78 

65 

61 

55 

As  the  non-bold  numbers  indicate,  TS  produced  more  temperatures  than  CV  in  8  out 
of  the  12  regions.  The  regions  varied  in  size  and  so  the  numbers  of  grid  points  varied 
from  region  to  region.  As  a  result,  weighted  averages  of  production  were  calculated.  The 
weighted  average  multiplies  the  CONUS  averages  by  six  and  the  Saudi  Arabian  averages 
by  two.  The  weighted  averages  for  TS  and  CV  were  49%  and  42%  respectively.  The 
non-weighted  production  averages  were  also  calculated.  The  percentages  were  the  same  - 
TS  produced  on  average  49  %  while  CV’s  average  production  rate  was  42  %.  The 
similarity  is  due  to  the  general  trend  of  larger  production  percentages  over  regions  with  a 
smaller  number  of  grid  boxes. 


53 


Even  though  TS  was  the  overall  production  winner,  TS  did  show  a  decline  in 
production  rates  from  summer  to  winter  for  all  regions.  The  TS  decline  was  expected  as 
the  snow  land  type  is  turned  off  for  TS  thus  reducing  the  temperature  production.  Even 
though  the  snow  land  type  is  turned  off,  TS  still  was  competitive  with  CV,  except  for 
Bosnian  and  Korean  Winters.  The  lack  of  TS  decline  into  winter  is  not  evident  in  Saudi 
Arabia,  which  is  a  region  naturally  hard-pressed  to  accumulate  snow.  Similarly,  the 
CONUS  showed  little  decline  in  production,  as  compared  to  Korea  with  the  largest 
decline,  as  the  majority  of  stations  in  CONUS  were  military  stations  concentrated  in  the 
Southern  CONUS,  again  where  snow  is  relatively  rare.  CV  did  not  have  a  downward 
trend  like  TS  from  summer  to  winter  except  for  perhaps  Korea.  However,  the  CV 
downward  trend  in  Korea  was  not  as  steep  as  TS’s  Korean  downward  trend. 

Looking  at  regional  aspects,  in  general  the  overseas  locations  had  more  temperatures 
from  both  algorithms.  More  data  in  potential  contingency  locations  is  definitely 
advantageous  for  US  military  operations.  Korea  had  both  the  largest  (summer,  TS)  and 
smallest  production  rates  (winter,  TS).  For  CV,  the  highest  production  was  in  the  Saudi 
Arabian  Fall  with  the  smallest  production  in  the  CONUS  Summer.  Saudi  Arabia  had  the 
highest  average  production  rate  for  aU  three  seasons.  This  is  a  bit  surprising  considering 
Saudi  Arabia’s  low  latitude,  which  results  in  the  smallest  satellite  coverage  of  the  three 
regions  (recall  Figure  3).  The  large  regional  variances  probably  reflect  the  differing  land 
types  from  region  to  region  and  the  algorithm’s  ability  to  correctly  determine  those  land 
types.  Saudi  Arabia  has  the  most  homogenous  land  type.  Since  the  algorithms  show  they 
could  identify  sandy  soils,  large  production  rates  were  understandable  in  this  region. 
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In  summary,  summer,  non-CONUS  sites,  and  TS  produced  more  temperatures,  but  a 
large  variance  of  production  rates  both  seasonally  and  regionally  was  evident.  On 
average,  TS  produced  7  %  more  data  than  CV. 

Data  Pair  Sample  Sizes 

A  data  pair  consists  of  the  ground  station  temperature  matched  with  a  TS  or  CV 
algorithm  produced  temperature.  Table  12  contains  the  number  of  data  pairs  by  region 
and  season. 


Table  12.  Number  of  Data  Pairs  for  the  Temperature  from  Sensor  Microwave  Imager 
(TS)  and  CalibrationA'^  alidation  (CV)  Algorithms 


SUM 

MER 

FALL 

WINTER 

TS 

CV 

TS 

CV 

TS 

CV 

CONUS 

1440 

651 

1157 

999 

1415 

1330 

BOSNIA 

26 

22 

26 

21 

222 

117 

KOREA 

949 

465 

731 

533 

434 

796 

SAUDI 

ARABIA 

546 

445 

390 

352 

152 

102 

Table  12  shows  TS  produced  more  data  pairs  (1 1  out  of  12  cases)  than  CV.  With  the 
production  data  in  mind,  one  would  indeed  expect  a  large  number  of  TS  data  pairs.  On 
the  other  hand,  with  the  same  reasoning,  one  would  expect  the  overseas  regions  to  have  a 
large  number  of  data  pairs.  However,  the  small  number  of  overseas  pairs  was  primarily 
due  to  a  less  dense  network  of  weather  stations  that  also  do  not  report  every  hour.  Many 
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foreign  stations  only  report  the  weather  every  three  hours.  Thus,  fewer  overseas  ground 
weather  stations  means  fewer  potential  matches. 

Although  not  evident  in  the  foreign  numbers,  the  CONUS  data  pair  numbers  dipped  in 
the  fall  and  rose  again  in  the  winter.  That  trend  was  not  entirely  evident  for  CV  because 
the  CV’s  summer  data  set  had  a  several  days  with  no  data.  The  fall  reduction  was  most 
likely  due  to  the  weather  station  reporting  code  change  (METAR)  which  occurred  at  this 
time  in  the  US  and  resulted  in  some  bad  and/or  lost  data  (Wall,  1997). 

The  Bosnian  samples  were  markedly  small.  The  small  sample  size  was  due  to  bad 
Date/Time  data  present  on  the  SSM/I  data  tape  sent  from  AFWA.  The  day  number  was 
incorrect  for  almost  all  the  days  used  in  the  Bosnian  region.  In  fact,  the  majority  of  the 
summer  days  for  Bosnia  had  the  same  day  number.  For  example,  the  21  Aug  96  file 
contained  the  Julian  day  for  5  Aug  96.  In  those  cases,  the  data  was  copied  into  the  correct 
file  assuming  the  Julian  day  matched  with  the  correct  satellite  observation  time.  This 
procedure  may  have  been  incorrect  as  wiU  be  pointed  out  in  the  next  two  sections. 

Bias 

Bias  (mean  error)  is  the  tendency  for  the  retrieved  temperature  to  be  too  high  (warm) 
or  too  low  (cold)  compared  to  reality  (station  data).  Although  not  as  important  as 
accuracy  to  AFWA  as  they  correct  for  bias,  bias  is  important  for  operational  users  of  CV 
surface  temperature  products.  Bias  was  calculated  for  each  region  by  season  and  is  shown 
below  in  Table  13. 
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Table  13.  Bias  (K)  for  the  Temperature  from  Sensor  Microwave  Imager  (TS)  and 
Calibration/Validation  (CV)  Algorithms. 

Bold  indicates  warm  bias.  Itahcs  indicates  CV  bias  >  TS  bias. 


SUM 

MER 

FALL 

WINTER 

TS 

CV 

TS 

CV 

TS 

CV 

CONUS 

-6.7 

-6.6 

-6.0 

-3.5 

-6.0 

-2.7 

BOSNIA 

-17.9 

-18.9 

-6.8 

19.7 

-1.2 

KOREA 

-10.6 

-7.4 

-8.6 

-5.1 

-5.3 

-2.8 

SAUDI 

ARABIA 

-6.6 

0.1 

-3.5 

2.3 

-2.8 

m 

The  boldface  indicates  warm  bias.  Table  13  shows  both  algorithms  had  a  cold  bias  in 
20  out  of  24  cases.  The  cold  bias  is  perhaps  to  due  to  the  effect  of  moisture  or  lack 
moisture  on  the  ground.  If  the  training  data  sets  were  made  over  land  that  was  more 
moist  than  normal,  an  average  dryer  land  condition  would  provide  less  microwave 
radiation  to  the  sensor,  resulting  in  a  smaller  algorithm  produced  temperature.  An 
analysis  of  soil  moisture  conditions  and/or  precipitation  days  was  not  accomplished.  As 
indicated  by  italics,  only  the  Bosnian- summer  contained  a  CV  bias  greater  than  the  TS 
bias.  In  all  other  cases,  the  CV  bias  was  smaller. 

For  the  majority  of  cases  for  both  algorithms,  the  bias  exhibited  a  seasonal  trend  of 
decreasing  from  summer  to  winter. 

From  a  regional  perspective,  Bosnia  clearly  had  the  largest  biases,  cold  or  warm.  The 
explanation  for  this  probably  was  due  to  the  erroneous  satellite  times  data.  The  small 
sample  sizes  probably  had  an  impact  on  the  large  biases.  However,  recall  the  Bosnian-fall 
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also  had  a  very  small  sample  size  and  its  biases  were  of  the  same  magnitude  as  those  of  the 
regions  with  large  data  sets.  It  is  interesting  to  note,  knowing  that  both  algorithms  use  the 
same  satellite  time  file,  that  the  Bosnian  Winter  produced  a  CV  bias  comparable  to  other 
locations,  but  the  worst  bias  of  the  whole  study  for  the  TS  algorithm  (19.7).  The  Bosnian 
data  is  certainly  suspect  due  to  the  errors  in  the  satellite  times. 

Another  regional  aspect  was  CV’s  warm  bias  in  Saudi  Arabia  for  all  three  seasons.  In 
contrast  to  the  previous  seasonal  trend,  the  CV  warm  bias  increased  from  summer  to 
winter.  An  explanation  for  the  warm  bias  is  possibly  due  to  the  different  land  types  in 
Saudi  Arabia  as  compared  to  the  Western  Plains  and  Dessert  of  the  US  where  NRL 
gathered  the  land  type  scheme  training  data.  To  produce  a  warm  bias,  the  Saudi  Arabian 
land  type  would  have  to  have  a  higher  emissivity  in  reality,  providing  more  microwave 
radiation  to  the  SSM/I  sensor,  and  in  turn,  providing  a  higher  algorithm  surface 
temperature  on  average. 

The  calculation  of  average  biases  for  season  and  region  was  planned  but  due  to  the 
trends  and  large  differences  between  the  algorithm  biases,  averages  were  not  calculated. 
Too  much  information  would  be  lost  in  the  generalization.  If  a  weather  observer  or 
forecaster  uses  a  TS  or  CV  surface  temperature  product  for  one  of  the  regions  and 
seasons  in  this  study,  they  should  consult  Table  13  when  no  other  guidance  is  available. 

In  summary,  both  algorithms  tended  to  produce  colder  temperatures  on  average,  with 
the  notable  exception  of  CV  in  Saudi  Arabia.  TS  temperatures  exhibited  the  largest 
departures  from  truth.  Most  biases  decreased  from  summer  to  winter. 
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Root  Mean  Square  Error  (RMSE)  -  Accuracy 

RMSE  will  answer  the  other  major  question  in  this  research,  “How  accurate  are  the  TS 
and  CV  algorithms?”  More  specifically,  did  the  algorithms  remain  below  the  AFWA 
eight-degree  cut-off  for  acceptable  accuracy?  Table  14  shows  RMSE  value  for  all  regions 
and  seasons.  Boldface  indicates  RMSE  greater  than  or  equal  to  8  K  (outside  of 
acceptable  accuracy  limits).  All  twelve  cases  were  included  in  the  analysis.  However,  the 
Bosnian  winter  was  suspect  due  the  large  difference  between  TS  and  CV. 


Table  14.  RMSE  (K)  for  the  Temperature  from  Sensor  Microwave  Imager  (TS)  and 
CalibrationA^alidation  (CV)  Algorithms.  Bold  indicates  RMSE  >  8  K. 


SUM 

MER 

FALL 

WINTER 

TS 

CV 

TS 

CV 

TS 

CV 

CONUS 

8.3 

7.8 

7.3 

5.6 

7.7 

6.2 

BOSNIA 

18.5 

19.4 

7.5 

5.3 

20.5 

5.2 

KOREA 

11.3 

8.0 

9.5 

6.2 

6.2 

6.2 

SAUDI 

ARABIA 

9.2 

■I 

■I 

8.0 

■I 

CV  was  more  accurate  than  TS  in  8  out  of  12  cases.  The  average  algorithm  RMSE 
accuracy  was  9.1  and  8.1  for  TS  and  CV,  respectively.  Thus  in  the  aggregate,  CV  was 
1  K  RMSE  more  accurate  than  TS.  There  was  no  clear  regional  or  seasonal  pattern  to 
algorithm  performance  as  all  the  differences  were  rather  close:  within  2  K  RMSE. 
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10  out  of  24  cases  failed  to  meet  the  accuracy  criteria  of  <  8  K  RMSE.  More 
importantly  for  AFWA,  TS  failed  to  met  the  criteria  in  6  out  of  12  cases.  Note,  TS  was 
close  (within  0.5  degrees)  to  the  8  K  RMSE  criteria  for  two  cases  (CONUS  winter  and 
Bosnian  fall)  and  CV  was  close  in  one  case  (CONUS  summer).  CV  failed  the  accuracy 
criteria  in  4  out  of  12  cases.  In  other  words,  CV  met  the  accuracy  (8/12  -6/12)  16  % 
more  often  than  TS.  Note,  if  the  accuracy  criteria  was  <  8  K  RMSE,  two  CV  scores 
(Korean  summer  and  Saudi  Arabian  fall)  would  be  changed  to  acceptable,  increasing  the 
acceptance  rate  by  50%. 

Seasonally,  in  summer  (excluding  Bosnia)  TS  was  less  accurate  than  CV,  and  in  fact 
never  met  the  criteria.  Summer  was  also  the  least  accurate  season  for  CV.  There  was  no 
consistent  trend  from  summer  to  winter  for  either  algorithm. 

In  short,  CV  was  more  accurate,  more  often  than  TS.  CV  met  the  accuracy  criteria 
16  %  more  often  than  TS  and  beat  TS  in  75%  of  all  the  cases.  On  average,  CV  was  1  K 
RMSE  more  accurate  than  TS. 

Summary 

In  terms  of  production,  TS  produced  7  %  more  surface  temperatures  than  CV.  On  the 
other  hand,  TS  was  less  accurate  than  CV.  On  average,  CV  was  1.0  K  RMSE  more 
accurate  than  TS  and  met  the  accuracy  criteria  16  %  more  often  than  TS.  Both 
algorithms  mostly  exhibited  a  cold  bias,  but  CV  displayed  a  warm  bias  in  Saudi  Arabia. 
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V.  Summary  and  Recommendation 


Summary 

With  a  large  comparative  study,  the  present  research  provided  information  on  AFWA’s 
two  main  questions  about  the  TS  and  CV  algorithms.  Concerning  amounts  of  data,  TS  on 
average  produced  7  %  more  surface  temperatures  than  CV.  The  production  percentages 
varied  greatly  by  region  and  season. 

The  accuracy  numbers  did  not  vary  as  much  by  region  and  season.  On  average,  CV  is 
16  %  more  accurate  in  meeting  the  AFWA  accuracy  criteria  of  RMSE  <  8  K.  CV  had  an 
average  1  K  RMSE  better  accuracy  than  TS.  CV  met  or  beat  TS  in  75  %  of  the  cases. 

It  is  interesting  to  note  the  bias  behavior  of  the  two  algorithms.  An  overall  cold  bias  is 
prevalent  but  warm  biases  did  exist  as  for  CV  in  Saudi  Arabia.  Biases  generally  declined 
in  winter. 

Data  Processing  Recommendations 

Several  or  all  of  the  eight  FORTRAN  program  types  could  be  combined  into  one 
program.  Certainly,  the  Reformat  Phase  (steps  B  &  C  and  H)  programs  could  be 
combined.  Combining  would  save  much  processing  time.  Similarly,  a  program  could  be 
written  to  automatically  combine  all  the  daily  match  files  into  the  single  file  with  all  the 
data  pairs  (last  part  of  Phase  4).  Doing  so  would  save  time  and  eliminate  the  manually 
intensive  cut  and  paste  operations  currently  needed  in  the  last  data  processing  step. 
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Research  Recommendations 


A  performance  comparison  should  be  made  of  over  the  various  land  types.  In  addition, 
the  effect  of  clouds  and  precipitation  on  algorithm  accuracy  could  be  investigated.  A 
correlation  study  could  be  done  to  match  up  weather  observations  of  surface  temperature, 
clouds  and  precipitation  with  algorithm  produced  temperatures. 

For  CV,  a  fifth  algorithm  is  available  which  is  a  combination  of  all  the  data  used  in  the 
regression  analysis.  The  fifth  algorithm  can  be  used  if  the  original  four  land  types  cannot 
be  determined  (McFarland,  1991).  The  effect  of  the  fifth  algorithm  on  CV  production  and 
accuracy  would  be  interesting.  Similarly,  contact  the  CV  contractor  should  be  done  to 
determine  if  any  new  CV  algorithms  for  other  land  types  are  available. 

Concerning  the  present  study’s  production  results,  CV  produced  temperatures  for 
approximately  90  %  of  all  the  grid  points  on  21  Oct  96.  No  other  day  for  either  algorithm 
had  such  a  high  rate.  A  cause  could  not  be  determined  and  is  left  for  further  study. 

Finally,  a  regression  analysis  could  be  done  on  this  study’s  data  pairs  to  derive  new 
coefficients  and  compare  them  to  the  current  coefficients  used  in  both  algorithms.  This 
study’s  database  is  by  far  larger  (>13,300  pairs)  than  the  studies  performed  at  AFWA 
(~1000  data  pairs)  (Kopp,  1997)  and  the  original  CV  validation  study  performed  in  1987 
(~1800  pairs)  (McFarland,  1991).  Separate  algorithms  could  be  generated  for  each 
location  and  season  analyzed  in  this  study. 
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Appendix  A  -  Summer  CONUS  Matches 


Appendix  B  -  Fall  CONUS  Matches 


Appendix  C  -  Winter  CONUS  Matches 


IVlOll  _  I  I  I  I,.  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  sAeg-g 


Appendix  D  -  Summer  Bosnian  Matches 
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Appendix  E  -  Fall  Bosnian  Matches 


Appendix  F  -  Winter  Bosnian  Matches 
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Appendix  G  -  Summer  Korean  Matches 


Appendix  H  -  Fall  Korean  Matches 


Appendix  I  -  Winter  Korean  Matches 
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Appendix  J  -  Summer  Saudi  Arabian  Matches 


Appendix  K  -  Fall  Saudi  Arabian  Matches 


Appendix  L  -  Winter  Saudi  Arabian  Matches 
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The  Air  Force  Weather  Agency  (AFWA)  has  two  operational  algorithms  that  derive  surface  temperatures  from 
microwave  observations  taken  by  the  Special  Sensor  Microwave  Imager  (SSM/I)  which  rides  aboard  space  platforms  of  the 
Defense  Meteorological  Satellite  Program  (DMSP).  The  algorithm  called  Temperature  from  Satellite  Microwave  Imager 
(TS)  is  used  to  analyze  global  cloud  coverage.  The  second  algorithm  is  fittingly  called  CalibrationA^alidation  (CV),  as  it 
was  the  algorithm  used  to  calibrate  and  validate  the  first  SSM/I  in  1987.  Multiple  linear  regression  defined  the  algorithms 
from  empirically  gathered  brightness  temperatures  and  simultaneous  surface  temperatures.  The  key  questions  were  how 
much  data  do  these  algorithms  produce  and  how  accurate  is  it. 

These  questions  were  answered  with  a  multi-seasonal  comparative  study  over  four  locations.  The  study  matched 
algorithm  outputs  to  conventional  weather  station  temperature  readings.  Over  13,300  data  pairs  were  generated  from  the 
1996  summer  and  fall  and  1996-1997  winter  for  the  Continental  United  States,  Bosnia,  Korea,  and  Saudi  Arabia. 

The  results  show  TS  produced  on  average  7  %  more  surface  temperatures  than  CV.  CV  met  AFWA’s  accuracy  criteria 
16  %  more  often  than  TS.  On  average,  CV  was  1.0  degree  Celsius  more  accurate  than  TS.  The  study  generated  bias  tables 
for  all  locations  and  seasons. 
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